ON REGULARITY OF A BOUNDARY POINT FOR HIGHER-ORDER 
PARABOLIC EQUATIONS: TOWARDS PETROVSKII-TYPE 
CRITERION BY BLOW-UP APPROACH 



V.A. GALAKTIONOV 

Abstract. The classic problem of regularity of boundary points for higher-order PDEs 
is concerned. For second-order elliptic and parabolic equations, this study was completed 
by Wiener's (1924) and Petrovskii's (1934) criteria, and was extended to more general 
equations including quasilinear ones. Since the 1960-70s, the main success was achieved 
for 2mth-order elliptic PDEs (e.g., by Kondrat'ev and Maz'ya), while the higher-order 
parabolic ones, with infinitely oscillatory kernels, were not studied in such great details. 

As a basic model, explaining typical difficulties of regularity issues, the ID bi-harmonic 
equation in a domain shrinking to the origin (0, 0) is concentrated upon: 

u t = -u xxxx in Qq = {\x\ < R(t), -1 < t < 0}, 

where R(t) > is a smooth function on [—1,0) and R(t) — * + as t — > - . The zero 
Dirichlet conditions on the lateral boundary of Qo and bounded initial data are posed: 

u = u x = at x — ±R(t), —1 < t < 0, and u(x, —1) = Uo(x). 

The boundary point (0, 0) is then regular (in Wiener's sense) if w(0, 0~) = for any data 
Mo, and is irregular otherwise. The proposed asymptotic blow-up approach shows that: 

(i) for the backward fundamental parabolae with R(t) = /(— i) 1 / 4 , the regularity of its 
vertex (0,0) depends on the constant I > 0: e.g., I = 4 is regular, while I = 5 is not; 

(ii) for R(t) = (-i) 1/4 ^(~ln(-t)) with tp{r) — > +oo as r — > +oo, regularity /irregularity 
of (0, 0) can be expressed in terms of an integral Petrovskii-like (Osgood-Dini) criterion. 
E.g., after a special "oscillatory cut-off' of the boundary, the function 

R(t) = 3-t2T(_t)i[ m |ln(-t)|] f 

belongs to the regular case, while any increase of the constant 3~i 2^" therein leads to 
the irregular one. The results are based on Hermitian spectral theory of the operator 
B* = -D { y i] - | yD y in L 2 p ,(R), where p*(y) = e~ a ^'\ a = constant G (0,3-2-f), 
together with typical ideas of boundary layers and blow-up matching analysis. Extensions 
to 2mth-order poly-harmonic equations in M. N and other PDEs are discussed, and a 
partial survey on regulariy/irregularity issues is presented. 



Dedicated to the memory of Professor f. G. Petrovskii 
Date: January 26, 2009. 

1991 Mathematics Subject Classification. 35K55, 35K40. 

Key words and phrases. Higher-order parabolic equations, boundary regularity, Petrovskii criterion, 
non self-adjoint operators, boundary layer, matching. Submitted to: NoDEA. 

1 



1. Introduction: Petrovskii's regularity criterion (1934) and extensions 

1.1. First discussion: regularity as a fundamental issue in potential theory. 

Well-posedness of initial-boundary value problems (IBVPs) for linear and nonlinear par- 
tial differential equations (PDEs) are key for general PDE theory and applications. The 
principal question, which was always in the focus of the key research in this area from the 
nineteenth century to our days, is to determine optimal, and as sharp as possible, condi- 
tions on the "shape" of the continuous boundary, under which the solution is continuous 
at the boundary points. This means that, for a given PDE with sufficiently smooth 
coefficients, a standard IBVP can be posed with a classical notion of solutions applied. 

It is impossible to mention all the cornerstones of such regularity PDE theory achieved 
in the twentieth century. In the nineteenth century, this study, for the Laplace equation 

(1.1) Au = f in O C t 2 or K 3 (0 G OVl is the point under scrutiny), 

began by Green (since 1828), Gauss (1940), Lord Kelvin (1847) and Dirichlet himself (in 
the 1850s), Weierstrass and Neumann (in the 1870s), Hilbert (1899), Schwarz, Poincare 
(since 1887); see detailed history of potential theory in Kellogg [Ml pp. 277-285]. It is 
also worth mentioning that classic in potential theory "Lyapunov's surfaces" were intro- 
duced about 1898 for general non-convex domains in R 3 (in the convex case, the problem 
was solved in 1870 by Karl Neumann). Lyapunov's study was known to be inspired by 
Poincare's earlier papers; a full proof for general Lyapunov surfaces was completed by 
V.A. Steklov about 1902. As a related issue, as was pointed out by V.G. Maz'ya [69] . 
Poincare seems was the first who in last years of the nineteenth century already knew 
and used the exterior cone condition for the regularity of boundary points. As is well 
known from his private letters, discussions with several Russian mathematicians, includ- 
ing e.g., his supervisor P.L. Chebysho\E| and Steklov, Lyapunov was constantly rather 
anxious studying new and amazing but not fully rigorously justified ideas and methods 
by Poincare. Possibly, this led him to develop the concept of "Lyapunov's surfaces" in 
potential theory about that time (as well as some of other fundamental concepts in his 
stability theory, where discussions with Poincare were also known to take place). 
For the second-order Laplace equation (jl.ip and the heat equation 

(1.2) Ut = Au in Q, (0,0) G dQ is characteristic, where dQ is tangent to {t = 0}, 

regularity theory was almost fully completed in the 1920-30s, and various ideas and key 
results on optimal regularity criteria created the amazing history, which is explained in a 
number of classic monographs. For convenience, we will begin our discussion with some 
historical aspects in next Section [31 where we hope to present some new and not that 
well-known regularity features and peculiarities for the sake of the attentive Reader. 



A.M. Lyapunov Master's Thesis "On Stability of Spheroidal Equilibrium Forms of a Rotating Fluid" , 
written under the supervision of Chebyshov, was defended in the S. -Petersburg University, 27th January, 
1885. 
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The present paper uses the asymptotic blow-up approach from react ion- diffusion theory 
to such optimal boundary regularity questions. It is shown that, for a number of higher- 
order PDEs with essentially oscillatory kernels, these can be treated by blow-up evolution 
via approaching a "singular" boundary point. Surely, this is not a principal novelty not 
only for the heat equation (11.21) (Petrovskii, 1934), but also in elliptic regularity theory, 
where, together with other techniques, rescaled linear operator pencils were shown to be 
key; we refer to Kozlov-Maz'ya-Rossmann's monographs [521 [53] and surveys in [27J [50j 
l68l 194"] for details and full history. In this connection, Kondrat'ev's seminal paper [19] 
for higher-order linear elliptic equation^ (see Kondrat'ev-Oleinik's very detailed survey 
of 1983 |50j for up to the 1980s and the most recent monograph [12] for a huge number of 
more modern extensions), thought not devoted to the regularity issues, represented a novel 
and involved use of spectral properties of operator pencils for describing the asymptotics of 
solutions near corner singularities. This analysis assumes a kind of an "elliptic evolution" 
approach for elliptic problems, which then is not well-posed in Hadamard's sense but can 
trace out the behaviour of necessary global orbits that approach the singularity point. 
For nonlinear elliptic equations, these ideas got various developments and diversions since 
the 1970s; see Veron's monographs [9~T| [92] and [10] for further history and key references. 

Classic Wiener's and operator pencil ideas were successfully developed for many higher- 
order elliptic equations with a typical representative: 

(1.3) (—A) m u = f in QcR N (m > 2), 

in the papers by Maz'ya with collaborators since the 1970s; see e.g., [71] and [671 1681 152| [53] 
as a guide. Though regularity point analysis for (II .3p revealed a lot of new phenomena 
and mathematical difficulties relative the classic ones for m = 1 (see comments later on), 
there are many papers on this subject, where essential progress was achieved, especially 
in the recent ten years. 

1.2. Higher-order parabolic PDEs: on some known regularity results. For 2mth- 
order parabolic equations, where the basic model is the canonical poly-harmonic equation 

(1.4) u t = —(—A) m u, with m > 2, 

regularity questions have been also addressed in the literature, however essentially less 
actively than for elliptic PDEs referred to above. Among results of classic solvability 
theory, which can be found in a number of monographs on linear and nonlinear equations, 
systematic approaches to the study of non-cylindrical domains with a characteristic vertex 
(0, 0), wherqj the tangent plane {t = 0} is horizontal in the {x, t}-space appeared already 

2 As explained in |49) . previous results on normal solvability of elliptic boundary value problems in 
domains with angular (in R 2 ) and conical (R N , N > 3) points were obtained by Fufaev (m = 1, N = 2, 
1960), Volkov (1962), Birman and Skvortsov (1962), Lopatinskri (1963). Earlier Nikol'skii's results [78] 
in the 1956-58 on boundary differential properties of functions in the Nikol'skii spaces H r v defined on 
regions in R 2 with angular points and hence criteria of such a regularity for the Laplace equation should 
be mentioned. 

3 Here and later on, we indicate only the results that are related to our blow-up setting to be explained 
shortly; these deep papers contain other involved conclusions. 
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in the 1960s. V.P. Mihailov papers in 1961 and 1963 [721 1731 [71] treated the case, where 
the parabolic boundary of the domain Qo lies below the characteristic plane t — and has 
a typical form of the corresponding "backward" paraboloid of the fundamental solution^ 
(in [731 EH) this shape is perturbed by factor that makes the paraboloid "sharper", i.e., 
becomes "more regular"; see further comments on that). Almost simultaneously with the 
fundamental "elliptic" paper [19 , Kondrat'ev in 1966 published another key "parabolic" 
paper [IS], which is less knownj. Here, Kondrat'ev dealt with the singularity issues 
for more general than (11.41) parabolic equations and various boundary conditions on a 
fixed fundamental paraboloid. Complicated asymptotics of solutions near such singular 
points were obtained, and, as a result, a solvability criterion was proved for existence 
of weak solutions in special weighted Sobolev-Slobodetskii-type spaces. In 1971, Feigin 
[25] extended the results to the case when, in any neighbourhood of the vertex (0,0), 
the boundary lies on both sides of the characteristic plane {t = 0} (see [HS] for Feigin's 
contribution to elliptic regularity theory in the 1970s). These deep results initiated further 
subsequent study, which nevertheless was not that exhaustive and sharp as in elliptic 
theory, and we will try to explain why (the answer seems easy: another level of difficulty) . 
The current standing of this regularity parabolic theory can be traced out further by the 
MathSciNet, but the author reports that getting a clear convincing view is not easy. 

1.3. Main goal: deriving sharp asymptotics of solutions near vertex depending 
on boundary geometry. As a consequence of the above discussion involving a number 
of strong previous results, obviously, optimal conditions of regularity /irregularity of the 
vertex (0, 0) for the poly-harmonic flow (II. 4p demand the following: 

Goal I: sharp asymptotic expansion of solutions near (0, 0) are necessary 
(1.6) for rather arbitrary backward characteristic (non-fundamental) paraboloids 
to get regularity ot the vertex via the geometry of the boundary nearby. 

In other words, we check whether it is possible to characterize a sharp asymptotic be- 
haviour of classic (strong) solutions near the characteristic vertex, when the paraboloid 
is not given by the backward variable of the fundamental solution, i.e., can take various 



^Q.v. more recent Mihailov's research on existence of boundary values of poly-harmonic functions for 
domains with smooth boundary; see references in [75] . 

5 It is useful to compare the total number of citations in the MathSciNet (October 2008): the elliptic 
paper [53] has the record 228 citations, while the parabolic one [JS] has 7 citations, i.e., 32 times less! 
However, the proposed scaling techniques [48l [49] leading to operator pencils are almost identical in 
both parabolic and elliptic papers (moreover, the parabolic pencil analysis |48j directly applies to some 
ultraelliptic PDEs). The asymptotic expansions derived by spectral theory near the vertex 0, 

(1.5) u(x,t) ~ J2(j k) r" 1 Aj '/ 2m (lnr) fc ipkj (Aj are eigenvalues of pencils) 

are similar for the elliptic (r = \x\) and parabolic (r = — t) cases (the ln fc r, k = 0, l,...,kj multipliers 
reflect finite algebraic multiplicities kj of Xj with associated generalized C°° eigenvectors {ipkj})- Note 
that, even for the heat equation, the required pencil spectral properties are not that easy to get, [J]. 
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forms. Loosely speaking, we are going to mimic a Petrovskii-type criterion for higher- 
order parabolic equations such as (jl.4p (we must admit that, literally, this goal is not 
and cannot be achieved literally, but a partial progress is indeed doable). Once such an 
asymptotic has been obtained, one can re-evaluate the solution near the vertex in any 
necessary weighted metric, eventually to check whether solution exists (is continuous at 
the vertex) in weaker sense. 

As usual, we find useful, in order to characterize our blow-up approach, to consider first 
the simplest ID bi-harmonic equation 



which even in the one x-dimension provides us with a number of not that mathematically 
pleasant surprises. It seems that the main difference (and the origin of extra difficulties) 
between the higher-order elliptic case (11. 3p and the parabolic one (I1.7P is that the later has 
always oscillatory kernel of changing sign infinitely many times, while for ( 11.31) . the kernel 
is less oscillatory and even remains positive in some dimensions even for m > 2. As is well 
known since Maz'ya-Nazarov's results for m = 2, N > 8 in 1986 [70], the changing sign 
kernels of elliptic operators produce new phenomena, where classic regularity techniques 
may fail and even standard conical points can be irregular (and singular, i.e., unbounded 
at the vertex). Note that Kondrat'ev's sharp estimates such as (I1.5P [49J can also indicate 
such a possibility, if sharp bounds on first pencil's eigenvalues are known. Therefore, a 
simple ID model (I1.7P becomes a key parabolic PDE with such an infinitely oscillatory 
kernel, so we concentrate on this equation in what follows. 

As a unified issue of the present analysis of higher-order parabolic and other PDEs 
(i.e., without any order-preserving or positive kernel features), we aim a sharp treatment 
of boundary regularity conditions by using spectral theory of related rescaled non self- 
adjoint operators or pencils of such operators. A full justification of some "approximate" 
regularity conditions are very difficult and lead to open problems. As an independent 
feature of the regularity questions, which underlines their difficult nature, as far as we 
know, even for the simplest higher-order elliptic or parabolic operators with real con- 
stant coefficients, there is a clear and natural difficulty for determining sufficiently sharp 
regularity conditions for linear operators with kernels of essentially changing sign. Of 
course, we exclude those operators for which the regularity is guaranteed by standard 
embedding results, or when the irregularity is induced by some symplectic-Hamiltonian 
properties, which do not allow any shrinking of the domain under consideration (e.g., the 
latter is true for Schrodinger operators^ iD t + (— A) m , where L^-norm is preserved). For 
such operators, Wiener's and other related capacity-like techniques, which assume kernel 
positivity, are naturally expected to fail, so new techniques are necessary. 

Thus, our goal is to develop an asymptotic method of blow-up regularity analysis, firstly, 
for the simplest ID bi-harmonic equation, and next extend to poly-harmonic ones (jl.4p . 
In general, together with other examples of PDEs, we follow the principle (as we have 





The author thanks I.V. Kamotski for a discussion that clarified this property. 
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seen, such kernels have been studied much less in the mathematical literature) 

(1.8) Goal II: to deal with kernels having infinitely many sign changes, 

that can occur not necessarily for parabolic PDEs only. We will show that a kind of 
Petrovskii-like integral "criterion" for regularity/irregularity of boundary points [EH [82] 
(1934-35) can be visualized. However, for operators with oscillatory kernels of changing 
sign, a standard deterministic analysis of such integral Osgood-Dini-type conditions via 
divergence/convergence (as for the second-order or other operators with more positively 
dominant kernels) is not available. Secondly, we discuss possible extensions of the method 
to other types of linear and nonlinear PDEs with the kernels like (11.81) . Meantime, we 
continue our survey of the history of crucial for us regularity issues and integral criteria. 

2. Elliptic equations as the origin of regularity theory: Zaremba 
(1911), Lebesgue (1913), and Wiener (1924) 

2.1. Laplace equation: early days of existence and nonexistence theory. The 

first criterion (a necessary and sufficient condition) of regularity of a boundary point 
was the Wiener famous one for the Dirichlet problem for the Laplace equation 

Au = f in a bounded domain Q C M 2 , G dfl [u = on dfl, for definiteness) , 

derived in 1924, |93j (the case of arbitrary N > 2 was also embraced). It is formulated in 
terms of a diverging series of capacities (measuring the thickness of the complement of Q 
near 0) of a discrete family of domains shrinking to the given point on the boundary; 
see extra details and further discussions in Courant-Hilbert [I5[ p. 306] and in Kellogg 
EE P- 330]. 

Concerning nonexistence for the Dirichlet problem (see [SSI E21 E3] for details), the 
first rigorous nonexistence example belonged to Stanislav Zaremba in 1911 [H2J p. 310] 
(for a plane domain whose boundary has isolated points, not connected; e.g., the surface 
comprised a sphere and its centre); see comments in $M p. 285]. The key non-solvability 
example was constructed by Lebesgue in 1913 [58, 59 jj who discovered striking examples 
of irregularity. E.g., he showed that, for a domain in M. 3 with the inner cone obtained 
by rotating the function x 2 — e -1 ^' 1 for %\ > small, about the Xi-axis (then the points 
with %2 > e~ l l xi do not belong to the domain, i.e., the domain has a sharp cusp turned 
inwards, a Lebesgue spine), the origin is irregular. Figuratively speaking, despite such 
a thin spine connection of the with the boundary, "effectively" the centre remains 
"separated" reminding earlier Zaremba's nonexistence construction. The same is true for 
%2 = , with any e > (0 is regular if it can be touched from outside by any cone 

given by rotation of x 2 = x\ for any k > fixed); see Petrovskii's text book [531 p. 325]. 

^A full Lebesgue account on the Dirichlet problem includes his important earlier paper [56j , where the 
variational method was employed in 2D to solving under very general conditions, and another key one 
[57], where "barrier" (in Poincare's terminology) techniques were employed; |59j is his final summarizing 
all key ideas and results. 
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Carleman's PhD Thesis |T3] of 1916 was also one of the first study of singularities of 
solutions of elliptic PDEs at non-regular boundary points; see also [oTJl Introd.]. 

Forty years later, new examples of nonexistence (somehow in lines with Zaremba- 
Lebesgue-Uryson's ideas) were developed in [471 p. 129] by Kondrat'ev. For the Dirichlet 
problem Am + u = in Q, he produced further examples of non-solvability for not simply 
connected domains for N = 2 (sufficient number of disjoint small discs are removed), 
and domains in M. 3 such as a ball with line segments or narrow cylinders removed to 
get sufficiently small interior diameter (one end stays on the surface for not violating 
connectedness). These constructions extend to arbitrary N and 2m, including the most 
delicate range N > 2m. By an m-capacity technique, [47] also studied the question of 
^-regularity (k < m) at a boundary point A meaning that any u G W^ has k — 1 
derivatives vanishing there asymptotically (up to a set of the measure o(a N ) in \x — A\ < a 
as a — > 0). 

2.2. Uryson (1924), Keldysh (1941), and Miusskaya Square. A nonexistence ex- 
ample similar to Lebesgue's one was independently constructed by P.S. Uryson (1898- 
1924) in the early 1920s; his paper |89j was published in 1925, i.e., after his tragic death. 
From electrostatic point of view, which is natural for such problems, the electricity "flows 
down" from the cusp, i.e., cannot be retained there. Besides taking a physics curse from 
his PhD adviser N.N. Luzin (Luzin's postgraduate course was completed by him in 1921), 
the physical intuition of Uryson was possibly created even before, when he, being in a 
gymnasium in Moscow, began to work at the Shanyavskii University in the Miusskaya 
Square, Moscow, under the supervision of the eminent Russian physicist P.P. Lasarev, 
and even published a paper on experimental research in X-radiation; see [6] (it seems 
less known that Uryson was the first, before Luzin, supervisor of A.N. Kolmogorov in his 
student's period). This emphasizes Uryson's outstanding physical motivation and under- 
standing in the very early days of his scientific carrier. Both the Shanyavskii University 
and the Physical Institute, founded in 1911 for the outstanding physicist P.N. Lebedev 
(who first measured "pressure of the sun light"), were erected in 1910-1912 by Russian 
architect I.I. Ivanov-Schitz by using funding from Moscow's merchants^. Later on, since 
1953, on the basis of Lebedev's Institute building, was developed the Keldysh Institute 
of Applied Mathematical founded and headed until his death in 1978 by M.V. Keldysh. 

Interestingly (and key for the whole comment), Keldysh made a significant and novel 
contribution to Dirichlet's problem theory in 1941 in his celebrated paper [41] (see also 
[4*0]). which treated irregular points, as well as the new notion of stability of points and 
of the problem inside the domain relative perturbations of the boundary (the notion of 
capacity is discussed in Ch. II, while Ch. Ill is devoted to Wiener's regularity criterion). 
In Ch. V, he proved that a stable boundary point is always regular, but the converse is 



8 This old building (as we will see, somehow related to boundary regularity) can be seen in "The World 
of Andrey Sakharov" (he first visited it in 1945), in http://people.bu.edu/gorelik/Publications 

9 The host institution for the author of this paper for almost twenty five years (as an excuse for paying 
possibly too much unreasonable attention to those historical events around it). 
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not true (Keldysh balls M were constructed having zero area of unstable points on dM 
and of positive harmonic measure; see [571 15] for a modern exhibition and extension). 

It is also worth mentioning that ten years later, in 1951, Keldysh published the related 
paper [12] on elliptic equations that are degenerated on the boundary, where particularly 
the novel idea that a part of the boundary may be free from any conditions was introduced 
the first. These ideas were later developed by Fichera (since 1956), Oleinik, Radkevich, 
Kohn, Nirenberg (in the 1960s), etc.; see the monograph [79] and [H] for a more recent 
standing. For a full collection, let us mention Keldysh's fundamental paper [13] (same 
1951 year!) on linear operator pencil theory with applications to elliptic PDEs. Therein, 
several novel ideas including difficult n-fold completeness questions of root functions, 
i.e., eigen and associated ones, for the Keldysh pencils with variable coefficients, Keldysh 
chains, etc. were introduced. Its sequel in the 1950s includes works by I.M. Glazman, 
M.G. Krein, V.B. Lidskii, and others; these Keldysh methods of [13] were later used by 
Browder, Agranovich, Krukovsky, Agmon, Egorov, Kondrat'ev, Schulze, etc.; see [18] for 
details. The paper |43j was later classified as one of the main achievements of pencil 
theory in the twentieth century (together with J.D. Tamarkin's results in his PhD Thesis 
in 1917, S.-Petersburgh's University [88]); see Markus' monograph [63] for the history and 
key results. 

2.3. On some recent achievements. Since we are not dealing with elliptic equations in 
what follows, it is worth mentioning here that a full extension of Wiener's regularity test 
(criterion) via the concepts of potential-theoretic Bessel (Riesz) capacities to 2mth-order 
strongly elliptic equations with real constant coefficients, 

(L(d)u = f, /GCg°(n), ueH™, 
(2- 1 ) ] L(d) = J2 a aP d a+ ^ a Ql3 = a Pw (-l) m L(0 > in R N \ {0}, 

^ \a\=\p\=m 

was completed in the case N = 2m in 2002 by Maz'ya0 (N > 2m also admits a treatment, 
though applied for a subclass of the so-called strongly elliptic equations; for iV < 2m — 1, 
the point is always regular by the classic Sobolev embeddings); see J67JEH] and references 
therein to earlier results and extensions. 

Concerning the related questions of sharp asymptotic behaviour near Lipschitz bound- 
ary at for problems (12. ip . in [51], a rather general case of the boundary given by 

tt = [ x = (x 1 , x n ) G R N : x n >ip(x')}, <p(0) = 0, 

is studied. The asymptotics as x — ► of a class of solutions are shown to depend on the 
Lipschitz constant of <f(y') in a neighbourhood of 0: 

k(p) = sup| y /| <p \V(p(y')\ (estimates are simplified if f Q ^-y^-dp < oo). 



The principal extension of Wiener's-like capacity regularity test to the nonlinear degenerate p- 
Laplacian operator — A p u = —V • (|Vu| p_2 Vm), p € (1, N], was also due to Maz'ya in 1970 [64]; later on, 
his sufficient capacity regularity condition turned out to be optimal for any p > 1. 
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Among others, let us also note an unusual Wiener's-type regularity and existence result 
for blow-up on the boundary solutions of the semilinear elliptic equation 

(2.2) Au-\u\ q ~ 1 u = in fi, u(x) -> +00 as x -> <9fi (g > 1); 

see [Ml E2] , where a long list of other related references is presented (equation (12. 2p will 
be cited again). 

3. Parabolic PDEs: Petrovskii's ^\og log criterion (1934) and around 

3.1. Parabolic regularity theory began with the heat equation: first results and 
definitions. From the beginning, for the second-order parabolic equations, the boundary 
point regularity analysis took a different direction (i.e., not of Wiener's type). This was 
Petrovskii [811 182] , who in 19340 completed the study of the regularity question for the 
ID and 2D heat equation in a non-cylindrical domain. For further use, we formulate his 
result in a blow-up manner, which in fact was already used by Petrovskii in 1934 |81j . 
This is about the question on a irregular or regular point (x,t) = (0,0) for nonnegative 
solutions of the IBVP 

{u t = u xx in Q = {\x\ < R(t), -1 < t < 0}, R(t) — >• 0+ as t — > 0", 
with bounded smooth data u{x, 0) = uq(x) > on [— R(— 1), R(— 1)]. 

Here the lateral boundary {x = ±i?(£), £ £ [—1,0)} is given by a function R(t) that 
is assumed to be positive and C 1 -smooth for all — 1 < t < and is allowed to have a 
singularity of R'(t) at t = 0~ only. We then study the value of u(x, t) at the end "blow-up" 
characteristic point (0, 0~), to which the domain Qq is said to "shrink" as t —>■ 0~. 

For the heat equation, the first existence of a classical solution result was obtained by 
Gevrey in 1913-14 [34] (see Petrovskii's references in [SU p. 55] and [H21 P- 425]), which 
assumed that the Holder exponent of R(t) is larger than |, i.e., in our setting, at t = 0~, 
all types of boundaries given by the functions 



(3.2) R(t) = (-*)" with any v > \ are regular (Gevrey, 1913-14). 



Note that, for 2mth-order parabolic equations, a similar result saying that, for 



(3.3) R(t) = (— £)2m, the problem is uniquely solvable 



(in Slobodetskii-Sobolev classes, i.e., not in the classic sense) was proved by Mihailov 
almost sixty years later. 

Regular point: as usual in potential theory, the point (x, t) = (0,0) is called regular (in 
Wiener's sense, see [67J; sometimes is also called boundary) if any value of the solution 
u(x,t) can be prescribed there by continuity as a standard boundary value on dQ$. In 



^Actually, this research was performed earlier: his most famous 1935 paper in Compositio 
Mathematica [55] was submitted for publication "(Eingegangen den 27. November 1933.)" [52J p. 419]. 
Before, the results were discussed, "Diskussion im wahrscheinlichkeits-theoretischen Seminar der Univer- 
sitat Moskau", which led to Kolmogorov's problem solved by Petrovskii in [82l p. 414]; see Abdulla [2j [3] 

for further details concerning this fundamental problem called nowadays Kolmogorov-Petrovskii's one. 
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particular, as a convenient and key for us evolution illustration, (0, 0) is regular if the 
continuity holds for any initial data Uq(x) in the following sense: 

(3.4) it = at the lateral boundary {\x\ = R(t), -1 < t < 0} w(0, 0") = 0. 

Irregular point: otherwise, the point is irregular (or inner) if the value u(0,0~) is not 
arbitrary and is given by the parabolic evolution as t — > CT (hence, formally and actually, 
(0,0) does not belong to the parabolic boundary of Qq, i.e., is inner in a natural sense). 

3.2. On some details of Petrovskii's analysis in 1934—35 and extensions. Reading 
last twenty five or so years various papers related to regularity issues for parabolic PDEs, 
the author found several and sometimes rather distinctive treatments of Petrovskii's re- 
sults. For instance, his 1934 paper [81] was not cited practically in no such papers, but 
the one [I]. There were also discrepancies in treating whether or not Petrovskii derived 
the integral Osgood-Dini-like conditions (surely, he did already in 1934). Therefore, we 
expect that it is worth mentioning here some aspects of Petrovskii's original derivation of 
his celebrated criterion. 

Thus, according to (13. ip . in both cases of regularity/irregularity, there occur interesting 
asymptotic problems of the behaviour of u(x,t) in Q as t — > 0~, which are of singular 
(finite-time blow-up) type. These problems were solved by Petrovskii [STJ in 1934-35, 
who derived a regularity criterion o for (13. 1 P by constructing rather tricky sub- and super- 
solutions ( "subparabolish" and " superprabolish" in the original Petrovskii terminology, 
[82] p. 386]) of the heat equation. 

For instance, to prove irregularity for the characteristic curve at (0,0") 

11 x 2 = 4(l+e)|t|ln|ln|t|| ", e > 0, 

he introduced the following explicit sub-solution "Barriere der Irregularitat" [811 p. 58], 
[82| p. 394] (see also comments in [76]): 

(3.5) " v(x,t)= ^-U, e^' + H^M " ( u * where i±a<l+ £ , 
and Ei < 1, k < 1 are some positive constants. Vice versa, to prove regularity for 

" x 2 =4|t|ln|ln|t|| ", (i.e., e = 0), 
the following explicit super-solution "Barriere der Regularitat" [8TJ p. 58] was used: 

(3.6) " v(x,t) = —±-r e^r + T? J = » (v t >v xx ). 

In [82] , the standard nowadays notation u(x, t) and u(x, t) for sub- and super-solutions 
were used. For more general boundaries given by 

(3.7) " x 2 < 4t\np(t) ", 

Petrovskii [821 P- 393] introduced the related function (p(t) determined from 

(3.8) " J^ = 61og^(t) ", 

*0 



Cf. quoting these necessary and sufficient conditions in Kondrat'ev [HI p. 450]. 
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which was used in constructing suitable super- and sub-solutions. Here, (13.81) is the origin 
of Petrovskii's criteria (13.111) given below. 

In his another paper [80] in 1934, Petrovskii also applied his novel method of sub- and 
super- solutions (i.e., before the classic Nagumo, or Miiller-Nagumo-Westphal, Lemmas 
from the 1940-50s; for the Dirichlet problem, barrier-type ideas in the methode de balayage, 
or sweeping out, were used by Poincare as early as in 1887 [M], with further development 
by Lebesgue in 1912 [57] and by Perron in 1923 for Am = [80], P- 432]) in the study of 
randomness that had a great influence on random processes theory, which he connected 
with linear elliptic equations. This role was first explained in Khinchin's monographs in 
1933 and 1936, [45] ; see below. Note that, in [SO], Petrovskii dealt with the Dirichlet 
problems ( "Dirichletschen Problems") for elliptic equations such as [80, p. 428] 

" M„£ + 2M xy £%- y + M yy + 2M.g + 2M y % = ", 

where ll M xx M yy > M xy " , and, surely, Petrovskii's typical double In In appeared as an 
its solution u u = lnln^J^=" on pp. 437-438. In three years time, sub-super-solution 

techniques were applied in the seminal KPP-paper, 1937 [46] . 

Petrovskii's 1934's paper [ST] is fully devoted to the optimal regularity and existence- 
uniqueness analysis for the heat equation, where he already introduced his 2vmln- 
criterion including the converging/diverging integral as in (13.111) for irregularity /regularity. 
For instance, he showed that the "5-log curve" and so on JSTJ p. 57] (cf. [S21 P- 404]) 

" x 2 = 4|t| [In | lnt|| + l|lnln| ln|t|| + In In | In |t||(l + e) In In In In | ln|t||] u. s. w. " 

is regular for e < and is irregular for any e > 0. In addition, for the 2D heat equation 

V>t u xx — u yy ^ 

it is mentioned [STJ p. 59] that the regularity /irregularity issues occur about the surface 

" x 2 + y 2 = —4(1 + e)£ln | In ". 

In particular, then [821 p. 397] the super-solution took the same form (13. 5p . where x 2 is 
replaced by x 2 + y 2 . In the later paper [82] in 1935, extensions to the heat equation in 
M fc x M + were mentioned. 

Hence, we attribute Petrovskii's criterion to 1934 [81], and not to 1935, [S2], as used had 
been done in most of other "parabolic" regularity papers, we have referred to here. Con- 
cerning related earlier papers on parabolic PDEs in non-cylindrical domains, Solonnikov 
(1965) [87] and Ivasishen (1969) [38] should be mentioned. 

Thus, using such novel barriers, Petrovskii JSU [82] established the following: 

(i) R(t) = 2v^t y/ln\\n(-t)\ (0, 0) is regular, and 

°" (n)R(t) = 2(l + e)V=t V ln l m (-*)l> £>° => (0, 0) is irregular. 

n 



More precisely, he showed that, for the curve expressed in terms of a positive function 
p(h) — » + as h — ► + ~ tj^t is about right) as follows: 



(3.10) R(t) = 2^i yj- lnp(-t), 

the sharp regularity criterion holds: 



(3.11) 



j p(h)y/\lnp(h)\ ^ < ( = ) + ^ ^ ^ Q _ ) ig irregular ( regular ). 



It is worth underlying again that both converging (irregularity) and diverging (regularity) 
integrals in (13.111) as Dini-Osgood-type regularity criteria already appeared in the first 
Petrovskii paper [EH p. 56] of 1934. It should be mentioned that, since Petrovskii used 
proposed him approach of constructing super- and sub-solutions, solving such partial 
differential inequalities led (as usual) unavoidably to a technical assumptions: for the 
super-solution (regularity of (0,0)), this is [H21 P- 392] 

" 2. *logp(*)-^ 0; " 

and, for the sub-solution (irregularity of (0,0)) [321 p. 397], 

"3. ilogp(t)->0 fur t^0; " and 

5. es ist ^ < gj^. 

These purely technical assumptions can be got rid of, pQ; our blow-up spectral-boundary 
layer approach is also assumptionless and rigorous for the heat equation, see Section 17.71 

The results of [81, 82] were key important for probability theory, where this result is 
expressed as follows: if Xt is a ID Wiener process and $(t) > is monotone such that 



OO 

(3.i2) jm e ~* 2 (t)/2 dt = 00j 

then, with unit probability, \x t \ < for all t ^> 1. On the contrary (this is about 
the criterion), if the integral converges, then with unit probability 3 {tk} — > oo such 
that \xt k \ > 3>(tfc); an optimal extension was done in [45] . Eventually, these Petrovskii's 
results led to the so-called the Law of Iterated Logarithms (LIL) [HI p. 392], which was 
discovered by Khinchin and Kolmogorov. Earlier, in 1924 Khinchin proved that, for a 
sequence {Xj}j>! of independent random variables with values ±1 with the probability 
|, there holds: 



(3.13) S n = X 1 + X 2 + ... + X n < y/(2 + e)nln\nn for all n > 1 (V e > 0) 

with the probability one. This sharp estimate improved earlier Hausdorf's inequality 
(1913) S n = o(n2 +£ ), Hardy-Littlewood's one (1914) S n < C\Jn Inn , and Steinhaus' 
inequality (1922) lim sup S n /\/2n In n < 1. As a final step, in 1929, Kolmogorov proved 
the LIL for any bounded independent random variables not assuming that the summands 
were identically distributed; see more details in [161 Ch. 7] and [35]. Clearly, (13.91) has 
probabilistic roots in (13.131) . 
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A Wiener-type criterion, already established by Khinchin [45] (an earlier result was in 
1933) in a probability representation, for the boundary regularity for problems such as 
( 13. ip and in WL N , was derived by Landis in 1969 [55J in terms of converging/diverging series 
of potentials of shrinking sets involved; see also [23] for a criterion via thermal capacity 
also obtained along the lines of derivation of Wiener's one for Laplace's equation. These 
results were not stated in the blow-up evolution and more "practical" Petrovskii's style 
(13.111) . though Wiener's-type capacity criteria serve for more general types of boundaries 
than according to Petrovskii's approach. 

Petrovskii's integral criterion of the Dini-Osgood type given in (13.111) is true in the 
iV-dimensional radial case with (see [3] [3] for the recent updating) 



a/| lnp(h)\ replaced by |lnp(ft,)| 2 . 

Several boundary regularity/irregularity results are now known for a number of quasilin- 
ear parabolic equations including degenerate porous medium operators. Though, some 
difficult questions remain open even for the second-order parabolic equations with order- 
preserving semigroups. We refer to [H [2j [31 [36J EH [76] as a guide to a full history and the 
extensions of these important results. 

As far as we know, (13.91) is the first clear appearance of the "magic" y/log log in PDE 
theory, currently associated with the "blow-up behaviour" of the domain Qo and cor- 
responding solutions. Concerning other classes of nonlinear PDEs generating blow-up 
y/log log in other settings, see references in [30] . 

Finally, as an introduction to difficult features of higher-order parabolic regularity to 
be studied, it is worth clearly stating that 



(3.14) for solutions of changing sign of the HE, Petrovskii's criterion ( 13. lip fails 



(obviously, positive subsolutions to prove irregularity are not applicable, while regularity 
can be still proved by comparison). In fact, a general criterion for arbitrary bounded 
changing sign data uq(x) in (13. ip cannot be derived in principle. The origin of this is the 
same for the HE (13. ip . the bi-harmonic one (11.71) . and others; see our analysis below. 

3.3. Example: a refined asymptotics for the porous medium equation. In con- 
nection with the blow-up log-log, it is worth mentioning the asymptotic result (seems, 
unique of this type) for the Dirichlet problem as in (13. ip [36[ § 2], which is now formu- 
lated for the radial porous medium equation in the pressure form: 

v t = Av m , m > 1, where the pressure u = v m ~ x solves: 

(3.15) 

u t = uAu + 7| Vw| 2 in Q 0} with R(t) = v^i, 7 = > 0. 

Then the log-log occurs in the asymptotic behaviour of nonnegative solutions as t — > 0~: 

(3.16) u(x,t) ~ 41n|1 * ( _ f)| ^0 uniformly on subsets - l| > ln|ln 1 ( _ f)| - 

The computations leading to (I3.16P by asymptotic-matching ideas (some will be involved 
in our further study) are not easy, and a complete proof was not supplied in [36] (regardless 
the comment at the end of page 4 in [36], it seems that a full and entirely rigorous 
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justification of this kind of matching with "floating" matching point can be extremely 
difficult and even seems illusive; this emphasizes a general complexity of such type of 
results even for second-order parabolic PDEs with the Maximum Principle). 

Remark: on boundary criterion for (13. 15j) - Obviously, the delicate asymptotics (13.161) 
allows one to detect the actual "criterion" of the boundary regularity for the equation 
( 13.151) (strangely, this was not addressed in [36]). We use the scaling invariance of ( 13. 15ft : 

(3.17) u = Aii, x = vAx, where A > is arbitrary. 

According to natural concepts of approximate similarity solutions (see e.g., [86, Ch. 6]), 
we next assume that 

(3.18) A = A(t) is not a constant, but a slow decaying function (relative to 

Then the equation for u(x,t) will include an extra very small asymptotic perturbation, 
which will not affect the asymptotics. Finally, according to (13. 16[) we set 

(3.19) A(t) = 41n |^ ( _ f) | , so that u ~ 1 

on the corresponding compact subsets as t — > CT. Thus, u(0, 0) ^ for the function 

(3.20) R(t) ~ 2 V /Z t ^/ln | ln( — ^t) | , 

and, in this irregularity condition, Petrovskii's magic vkln again mysteriously occurs. 
The constant "2" here seems not relevant as in (I3.9p . since (13. 15j) is nonlinear and admits 
also the symmetry u \— > Au, t i— > 4 for any A > 0, so that the domain behaviour 

with R(t) ~ \/^-t yJliL | ln(— 1)\ is irregular for (0, 0). 

* * * 

Thus, since the 1930s, Petrovskii's regularity \/log log-factor entered parabolic theory 
and generated new types of asymptotic blow-up problems, which have been solved for a 
wide class of parabolic equations with variable coefficients as well as for some quasilinear 
ones. Nevertheless, such asymptotic problems are very delicate and some of them of 
Petrovskii's type remain open even in the second-order case. 

For higher-order poly-harmonic operators, the situation becomes much more difficult, 
since the Maximum and Comparison Principles and order-preserving properties fail, so 
classic barrier techniques associated with sub-, super-, and other solutions are not appli- 
cable in principle. This precisely reflects the area of the present research: no Maximum 
Principle tools, so that the boundary regularity problem in Petrovskii's setting falls into 
the scope of a blow-up asymptotic behaviour study. 

4. Towards Petrovskii's criterion for the bi-harmonic flow: asymptotic 
blow-up setting, method, and layout of the paper 

4.1. The basic Dirichlet problem under consideration. As a basic model, using the 
"minimal" extension of the heat equation (13. ip by increasing the order of the parabolic 
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operator by two, we consider the bi-harmonic equation in the same shrinking domain as 
in (13. ip with the zero Dirichlet boundary conditions on the lateral boundary 8Qq: 



(4.1) 



u t = -u xxxx in Q = {\x\ < R(t), -1 < t < 0}, 
u = u x = at x — ±R(t), -1 < t < 0, 
u{x, 0) = u (x) on [-R(-1),R(-1)), 



where Uq(x) is a bounded and smooth function, uo(±R(— 1)) = 0. 

First of all, obviously, in modern PDE theory, there are many very strong regularity 
results, proved for wide classes of linear and nonlinear equations including, of course, 
the simplest among others bi-harmonic one (14. ip . We have mentioned some key results 
and papers before and will not indent to compete with those nice results and present no 
further references to this classic literature. The only issue we propose here is as follows: 
to show how to derive 

(4.2) a sharp boundary regularity criterion for (14. ID as a blow-up problem. 

As we have seen, these blow-up aspects of the regularity problem were already clearly 
revealed by Petrovskii in 1934. Though, for the second-order case, there are other equiva- 
lent approaches based on the positivity of the kernel and various order-preserving features, 
which fail for oscillatory kernels. 

To be more precise, this blow-up study for (14.11) is performed for: 

(4.3) a generic class of solutions exhibiting a "centre subspace" behaviour. 

In other words, as our Hermitian spectral theory of higher-order non self-adjoint linear 
operators in Section [5] shows, one can expect that, in addition to (14. 3p . 



It is principal for us that, to confirm the optimal character of the class of generic 
blow-up asymptotics, we will check (Section 17.71) that, by elementary spectral calculus, 



Hermite spectral theory yields precisely Petrovskii's criterion in (13.111) . 

Note that the problem: (14.31) or (14. 4p . also exists for the heat equation (13.11) (cf. (13. 14|) ) . 
but here, by the Maximum Principle, positive solutions always belong to the generic class 
(14. 3p . For the bi-harmonic PDE (14.11) . in a general setting, distinguishing generic patterns 
(14. 3p from those in (14.41) is difficult and seems even impossible (or makes no sense). 

4.2. Slow growing factor <p{r). Thus, similar to (13. 9p . we need to assume that 



(4.4) 



3 an infinite-dimensional set of "stable manifold" patterns, 
for which the conditions of regularity of (0, 0) are all different. 



(4.5) 



the same asymptotic approach to the classic problem (13. ip by classic 



(4.6) R{t) = (-*)* <p(r), where r = - ln(-i) +oo as t ->• 0" 
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Here, the scaling main factor (— t) 1 / 4 naturally comes from the bi-harmonic kernel variables 
(see (14.111) and (15.31) ). and <£>(r) > is an unknown slow growing function satisfying^ 

(4.7) <p(r) -> +oo, (p'(r) -> 0, and ^ -> as r -> +oo. 

Moreover, as a sharper characterization of the above class of slow growing functions, we 
use the following criterion: 

(4.8) 00 as r ^ +0 ° (^)^°)- 

This is a typical condition in blow-up analysis distinguishing classes of exponential (the 
limit is 0), power-like (a constant ^ 0), and slow-growing functions. See (86], pp. 390-400], 
where in Lemma 1 on p. 400, extra properties of slow-growing functions (14.81) are proved. 
For instance, one can derive the following comparison of such <p(r) with any power: 

(4.9) for any a > 0, (p(r) < r a and (p'(r) < r Q_1 for r > 1. 

Such estimates are useful in evaluating perturbation terms in the rescaled equations. 

Thus, the monotone positive function <p(t) in (14. 6p is assumed to determine a sharp 
behaviour of the boundary of Qq near the shrinking point (0, 0) to guarantee its regularity. 
In Petrovskii's criterion (13.91) . the almost optimal function, satisfying (14.71) . (14.81) . is 

(4.10) V ? i( r ) ~ 2Vlnr as r — > +oo, 

a dependence we have to compare our final results for the bi-harmonic equation with. 

4.3. First kernel scaling and layout. By (14. 6p . we perform the similarity scaling 

(4.11) u(x,t) = v(y,r), where y = j-^. 

Then the rescaled function v (y, r) now solves the rescaled equation 

{V T = B*V = -Vyyyy ~ \ yVy Ul Q Q = < If (t) , T > 0}, 

v = Vy = at y = ±<p(t), t > 0, 
v(0,y)=v (y) = u (y) on [-^(-1), R(-l)}. 

In view of the divergence (14.71) . it follows that our final analysis will essentially depend 
on the spectral properties of the linear operator B* on the whole line R. We reflect 
this Hermitian spectral theory in the next Section [SJ with application to the regularity 
criterion in Section [3 In particular, this differs our analysis from Kondrat'ev's classic one 
08]. 

However, in Section [6], we begin the study of the regularity of the vertex (0, 0) of the 
bi-harmonic backward fundamental parabolae: 

(4.13) R(t) = l(-t)*, i.e., <p(r) =1 = const. > 0. 

13 As we have mentioned, for the simpler case y(r), <p'{r) — > as r — > oo, the regularity was already 
proved by Mihailov in 1963 [73l [74]; in a certain sense, this extended the Gevrey-like result (|3.2|) for 
m = 1; see (|3.3p . 
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Then the problem (14.121) is considered on the fixed interval = {\y\ < /}, so that the 
final conclusion entirely depends on spectral properties of B* in Ii with Dirichlet bound- 
ary conditions. Of course, then the spectral problem for B* (not a pencil) becomes a 
very particular case of the general setting developed in [IB], but nevertheless, the clear 
conclusion on regilarity/irregularity becomes rather involved, where numerics are neces- 
sary to fix final details. In addition, as we pointed out, in more general setting for the 
fundamental backward paraboloids in M. N , the existence, uniqueness, and regularity of 
solutions in Sobolev spaces was proved in a number of papers such as [72j E3 EH [25], etc. 
Note that in [73j P- 45], the zero boundary data was understood in the mean sense (i.e., in 
the L 2 -sense along a sequence smooth internal contours "converging" to the boundary). 
Nevertheless, we have to stress attention to this simple case in order to reveal the exact 
transition between regularity and irregularity in the classic sense in the critical case f )4.13p . 
It seems that such a border case in between was not addressed before sharply enough. 

Thus, our conclusion is rather disappointing: unlike the classic heat equation (13.11) . for 
the bi-harmonic equation (14.11) . the vertexes of the fundamental parabolae ( 14. 13[) are not 
necessarily regula'^j. For instance, on the basis of careful numerics, we show that: 



R(t) 


= 4(-t)i = 


=> (0, 0) is regular, while 


R(t) 


= h{-t)l = 


=^> (0, 0) is irregular. 



Moreover, in the latter case, the vertex (0,0) is singular, i.e., 

(4.15) in general, u(x,t) is unbounded as (x,t) — > (0, 0~). 

We also claim that (by a continuity argument) there exists some l\ G (4, 5) such that 
(0,0) is irregular, but not singular, i.e., the limit u(0, 0~) exists. 

These regularity variations for constant functions v 5 ( r ) = I inevitably impose a special 
restriction for the study in Section [7J unbounded <^(t)'s satisfying (14.71) . where we derive 
a Petrovskii-like "criterion" of regularity. Namely, the irregularity condition is obtained 
in terms of an Osgood-Dini-type integral condition, while, in an accordance with (14. 141) . 
for the regularity, a special "oscillatory cut-off" of the lateral boundary must be applied. 

In Section [HI we discuss the extensions of the boundary regularity analysis to 2mth- 
order poly-harmonic operators 

(4.16) u t = (-l) m+1 D 2 x m u in Q (m>2), 

with the zero Dirichlet conditions, as well as for the same iV-dimensional problem for (11.41) 
in Qo, where the shape of the shrinking domain Qq is also of importance for the boundary 
regularity. For instance, one can consider the radial Qo given in (14.11) . where \x\ stands 
for the radial spatial variable. We also discuss regularity conditions for the third-order 

It is worth comparing this with the following well-known (since 1986) negative conclusion in elliptic 
theory: for fourth-order elliptic equations (with N > 8; for 5 < N < 7, the biharmonic capacity settles 
the regularity result in Wiener's sense [65] : N < 4 = 2m is done by Sobolev embedding), the vertex of a 
cone can be irregular-singular (the solution unbounded at the vertex), [70] . 
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linear dispersion equations, for a quasilinear fourth-order porous-medium-type equation 
(the PME-4), and for the linear wave (beam) equation of the fourth order. 

The present paper aims to give a first insight into principal difficulties of a sharp study 
of boundary regularity for higher-order parabolic and other evolution PDEs. Such a study 
inevitably generates a number of difficult mathematical problems, which, partially, do not 
appear in the second-order case, or can be avoided by positivity kernel properties. We 
must admit that some of them are not and even cannot be completely rigorously solved. 
The main problem of concern is that distinguishing generic patterns in (14. 3[) from non- 
generic ones in (I4.4p is not possible in general. Here, delicate boundary layer theory is 
essentially involved, where we ought to use some accurate numerical calculations by the 
MatLab whenever necessary to avoid huge technical digressions (using those seems also 
inevitable for truly 2mth-order parabolic equations with large enough m > 2). Regard- 
less a certain rigorous incompleteness of the mathematical analysis in some steps (which 
seems to be inevitable in general), we decide to demonstrate the whole machinery of the 
asymptotic blow-up methods in parabolic boundary regularity theory, and hope that this 
will help to attract some attentive Readers to improve the results when necessary. It 
is worth mentioning that, as can be seen in many directions of modern PDE theory of 
linear and nonlinear equations, the transition to higher-order models is accompanying by 
a dramatic increase of the complexity of the methods involved to achieve in a rigorous 
manner the necessary (basic or not) results. This increase of complexity sometimes mea- 
sures in orders, and often desired rigorous proofs can be illusive in a sufficient generality. 
A suitable restriction of the generality to achieve rigorous conclusions, though being an 
important research step, can lead to an extended and often unjustified number of artifi- 
cial hypothesis and even yield non-constructive assumptions (i.e., those that cannot be 
checked in a reasonable finite time or at all). In what follows, we prefer to avoid stating 
such theorems with non-constructive hypothesis, though this is indeed plausible in a few 
placeJ^l. In this case, we believe that the exchange of the ideas and methods, even in 
the case of a certain lack of a completely rigorous justification, can be key for further 
improvements and developing more consistent mathematical PDE theory. Of course, this 
imposes certain restrictions on the paper style, but the author hopes that the interested 
Reader will easily distinguish the rigorous and non-rigorous arguments to be used. 

5. Fundamental solution and Hermitian spectral theory for {B, B*} 

For convenience of the further study, we consider the 2mth-order poly-harmonic equa- 
tion (14.161) . and will describe the necessary spectral properties of the linear differential 
operator (the analogy of that in (14.121) for any m > 2) 



The main goal of a mathematician is not proving a theorem, but an effective investigation of 
the problem...", A.N. Kolmogorov, 1980s (the author apologizes for a non-literal translation from the 





Russian). 



IS 



and of its adjoint B in the standard L 2 -metric given by 

(5.2) B = (-l) m + 1 Dj m + ^ yD y + ± I (I denotes the identity). 

Both operators are not symmetric and do not admit a self-adjoint extension. We will 
follow [T7] in presenting necessary spectral theory. 

5.1. The fundamental solution and its sharp estimates. We begin with determin- 
ing the spectrum and the eigenfunctions of the adjoint operator B, which appears in 
constructing the fundamental solution of (I4.16P that takes the standard self-similar form 

(5.3) b(x,t)=t-&F(y), j/ = ^b. 

Substituting b(x, t) into (I4.16P yields that the radially symmetric profile F(y) is the unique 
even integrable solution of the linear ODE 

(5.4) BF = in R, J R F = 1, 

so it is a null eigenfunction of B. Taking the Fourier transform leads to 
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(5.5) F(y) = ceo f e s ™ cos(sy) ds, 



o 



where ojq > is the normalization constant, and, more precisely [20J, 



(5.6) F (v) = 7k J e Vs\y\ J -i(s|y|)ds in M, 

o 

where J v denotes Bessel's function. The rescaled kernel F(y) then satisfies a standard 
pointwise exponential estimate [19] 

(5.7) \F(y)\ < De- d °M a in R, where a = 

and D and d are some positive constants (d to be specified below). Such optimal 
exponential estimates of the fundamental solutions of higher-order parabolic equations are 
well-known and were first obtained by Evgrafov-Postnikov (1970) and Tintarev (1982); 
see Barbatis [HI [9] for key references and results. 

However, as a crucial issue for the further boundary point regularity study, we will 
need a sharper, than given by (15. 7J , asymptotic behaviour of the rescaled kernel F(y) as 
y — ► +oo. To get that, we re-write the equation (15.41) on integration once as 

(5.8) (_ 1 )m+l_ F (2m-l) + J_ y p = o j n R _ 

Using standard classic WKBJ-type asymptotics, we substitute into (15.81) the function 

(5.9) F(y) ~ y' So e ay " as y -> +oo, 
exhibiting two scales. This gives the algebraic equation for a, 

(5.10) (-1)1^ = ^, and 5 = ^- 1 >0. 

Note that the slow algebraically decaying factor y~ So in (15. 9p is available for any m > 2 
and is absent for m — 1 for the pure exponential positive Gaussian profile F; see below. 
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By construction, one needs to get the root a of (I5.10p with the maximal Re a < 0. This 
yields (see e.g., E] and [521 P- 141]) 

(5-H) a=g-i[cos(^ T )+isin(^ T )] ^-rfo + i&o (do = sin (5^) > 0). 
Finally, this gives the following double-scale asymptotic of the kernel: 

(5.12) F(y) = y- 5 °e- doya [C 1 sm(boy a ) + C 2 cos(b y a )]+ ... as y -> +00, 

where C\p. are real constants, \C\\ + | C 2 | 7^ 0. In (15. 12ft . we present the first two leading 
terms from the m-dimensional bundle of exponentially decaying asymptotics. 
In particular, for the bi-harmonic equation ( 14. lft . we have 

(5.13) m = 2 : a — |, <i = 3-2~^", 6 = 35-2~^ L , and 5 = |. 

5.2. The discrete real spectrum and eigenfunctions of B. We describe the spec- 
trum cr(B) of B in the space L 2 (R,) with the exponential weight 

(5.14) p(y) = e*M a > in R (0 = 5^), 

where a < 2do is a positive constant. Denoting by (•, -) p and || • || p the corresponding inner 
product and the induced norm respectively, we introduce a standard Hilbert (a weighted 
Sobolev) space of functions H 2m (M.) with the inner product and the norm 

2m 2m 

(v,w) p = fp(y) £ D k v{y)D k w{y)dy, \\vf p = j p{y) £ \D*v{y)\ 2 dy. 

K fc=0 E k=0 

Then F^ W (E) C L 2 p (R) C L 2 (M), and B is a bounded linear operator from H 2m (R) to 
L 2 (W). With these definitions, the spectral properties of the operator B are given by: 

Lemma 5.1. (i) The spectrum o/B comprises real simple eigenvalues only, 

(5.15) a (B) = {\ k = -±, k = 0,1,2,...}. 

(ii) The eigenfunctions ipk{y) ore given by 

(5-16) My) = ^D k y F(y) 

and form a complete set in L 2 (R) and in L 2 (R). 

(iii) The resolvent (B — A/) -1 for A (jL f(B) zs a compact integral operator in L 2 (R). 
The operators B (B*) have zero Morse index (no eigenvalues have positive real part). 

5.3. The polynomial eigenfunctions of the operator B*. We now consider the op- 
erator (15.11) in the weighted space L 2 *(M.) ((•, -) p * and || • || p . are the inner product and the 
norm) with the exponentially decaying weight function 

(5-17) p*(y) = ^ = e-«\y\ a > 0, 
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and ascribe to B* the domain i?5™(R), which is dense in iA(R). Then B* : H 2 p T 
L 2 , (R) is a bounded linear operator. Hence, B is adjoint to B* in the usual sense: denoting 
by (•, •) the inner product on L 2 (R), we have 

(5.18) (Bv,w) = (v,B*w) for any v G H 2m (R), w G #J m (R). 

The eigenfunctions of B* take a particularly simple polynomial form and are as follows: 
Lemma 5.2. (i) a(B*) = a(B). 

(ii) The eigenfunctions ipl(y) of'B* are polynomials in y of order k given by 

[k/2m] 

(5.19) r k (y) = ^ E 4P^ m V, A; = 0,1,2,..., 

3=0 

and form a complete subset in L 2 * (R) . 

(iii) B* has compact resolvent (B* — A/) -1 in L 2 *(R) for A ^ er(B*). 

With the definition (15. 19ft of the adjoint basis, by integrating by parts, the orthonor- 
mality condition holds (5m is the Kronecker delta): 

(5.20) 0>k,^) = hi for any k, I > 0. 

For m = 2 (this case will be treated in greater detail), the first eigenfunctions are 

V>o(y) = i, ^i(y) = j/, V>2(y) = 75 2/ 2 , V>3(y) = t^i/ 3 , 

(5 21) 

^(y) = ^3 (y 4 + 24), ^ B (y) = (y 5 + 120y), ^(y) = ^ (v* + 360^), 
etc., with the corresponding eigenvalues 0, —h, — ~, — |, —1, — |, — |. 

6. The vertex of fundamental parabolae can be regular or irregular 

Thus, consider the backward fundamental parabolae given by ( I4.13p . We continue 
denote by B* the corresponding linear operator 

(6.1) B*v = -v^-\yv' in I h v = v' = at y = ±Z, 

with the standard definition of the domain, etc. Here, B* is a regular ordinary differential 
operator with bounded smooth coefficients and a discrete spectrum cr(B*) [77]. Note that 
B* is not symmetric and does not admit a self-adjoint extension, so that er(B*) is not 
necessarily real (though we will observe a lot of real eigenvalues to be explained by a 
branching phenomenon). By A&(Z), we denote /-branches of eigenvalues of B* depending 
on the length I > 0. In view of Section [5j the limit I — > +oo is of a particular interest. 
The regularity criterion of the vertex is now easy: 



(6.2) vertex (0, 0) of parabolae (131511 is regular iff cx(B*) C C_ = C D {Re A < 0}. 



Then the solution v(y,r) (and any u(x,t)) gets exponentially small: as r — > +oo, 

(6.3) u(y, r) ~ e A ° T -> => |m(x, t)| ~ 0(r RoA °) -> as t -> 0", 

where Ao(Z) is the eigenvalue of B* with the maximal negative Re (in fact, Xo(l) turns out 
to be always real). Let us now briefly describe main steps of the current study. 
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6.1. The vertex (0,0) is regular for not that large I < 3.9779... . This conclusion is 
simple: as a standard practice, by multiplying in L 2 (Ii) the eigenfunction equation 

(6.4) b*#* = -(^*)( 4 ) - 1 y (m*y = \m* 

by the complex conjugate and summing up with the result of multiplying the complex 
conjugate equation by on integration by parts, one obtains 

(6.5) (A + A) / |^*| 2 = -2/ \Wf + \ J \^*\ 2 . 
Using the Poincare inequality: 

(6.6) J\V'\ 2 > Ao(0/l*| 2 , *e#oW; Ao(0 = fA (1), A (l) = 31.3..., 
where A (l) is the first eigenvalue of Dy > in Hq(Ii), we obtain from (16. 5ft : 

(6.7) (A + A) / |^| 2 < - -\]J |^| 2 . 

Therefore, the regularity of the vertex (0, 0) is guaranteed for (this is not sharp but close): 

(6.8) ^W_J>o l<k = [8Ao(l)]' = 3.9779.... 

6.2. First eigenvalues of B* are real. Here and later on in a few special technically 
difficult cases, we are going to rely on a clear numerical evidence for our conclusions. To 
this end, we use the MatLab solver bvp4c with the enhanced accuracy and tolerances 

(6.9) from 10~ 10 to the minimal admitted 10 -13 . 

This will guarantee that real eigenvalues, which can be very small, are computed correctly. 

In Figure (TJ we demonstrate first eigenfunctions of the operator B* in (16. ip for 1 = 1, 
where \I/q, and ^ are even functions, while and ^ are odd ones. Note that, by 
obvious reasons, the first eigenvalue 

A (l) = -31.16... 

is rather close to the eigenvalue — A (l) = —31.3... in (16.61) of the self-adjoint counterpart 
—D y 4 \ i.e., the non-symmetric perturbation — \^yD y is negligible. Actually, the "Sturmian 
geometric structure" of eigenfunctions in Figure [TJ each ^%{y) has precisely k zeros in Ji, 
which is true for —D y ^ [22], clearly remains valid for B*. Moreover, we have observed that 
the first eigenfunctions o/B* and —D^ for I = 1 coincide within the accuracy 0.5 x 10~ 2 . 
In fact, this suggests to get Xk(l) by branching at fi = + from eigenfunctions of —D y ^ < 
by using the following family for /x 6 [0, J (a homotopy path) of operators in Hq(I{): 

(6.10) B* = —Dy 4 ^ — fxyD y = B* for p = \. 

For /i = 0, the self-adjoint Bq = —D y ^ has the desired real spectrum and complete-closed 
set of eigenfunctions; see [39] and [90] for classic perturbation/branching theory. Later on, 
we present another branching explanation of the origin of real eigenvalues/eigenf unctions 
of the non self-adjoint operator B* by using results from Section [5] reflected / = oo. 

22 



Five eigenfunctions of B : A, =-31 , ^ =-237, 3, \ 3 =-2497, \=-5573 

¥*(y) 




FIGURE 1. Five first real eigenfunctions of the operator (|6.ip for I = 1. 



Shooting first zero of X (I): X (4.078.. ,)=0 



*(y) 




<F(y) 



Shooting second zero of X (I): X (7.25...)=0 




(a) first zero Zi = 4.08... (b) second zero h — 7.25... 

FIGURE 2. Shooting zeros of Ao(0 : nrs t zero ( a ) an d the second one (b). 

6.3. Ao(Z) changes sign: numerical evidence. It is clear that, for large I, the non- 
symmetric perturbation —^yD y becomes essential in B*, and Sturm's zero property in 
Figure [1] does not remain true. In Figure [2(a), (b), we present the shooting results of 
numerical calculating first three roots of the branch of the first eigenvalue {Ao(/)}: 

(6.11) A o (0 = at first three roots h = 4.08... , l 2 = 7.25... , and l 3 = 10 

Note that li is about 3% close to the bound I* = 3.9779... derived in (|6.8p . Further results 
of numerical experiments are shown in Table 1 fixing these 3-4 roots of Xq(1). 

Thus, taking into account these numerical results, in view of ( 16. 2p . we conclude that 



(6.12) 



(0,0) is regular for I e (0, Zi), is irregular (singular) for / e [Zi,^], etc 



23 



Table 1. Some values of Xo(l) 



I 


U V / 


1 
1 


— 01.10... 


Z 


— 1.00... 


Q 

a 


n oka 7 


A 
-i 


— U.UUolOZ... 


ZL (17^ 


— u.uuuzou. . . 


A D77^ ~ / 


n nnnni 1 1 

U.UUUUllO... 




u.uuuzo. . . 


A 1 
4t. 1 


n DD99 


zL 9 


n ni 1 

U.U 1 J 


r. 
■J 


U.U^toO... 


c 




u.uioo.. . 


A 
U 


n D4R 


7.25 « Z 2 


0.00167... 


7.5 


-0.0097... 


8 


-0.027... 


9 


-0.018... 


10«Z 3 


-0.00084... 


11 


0.00397... 


12 


0.00172... 


13 « Z 4 


-0.00055... 



In addition, since Ao(^i) = 0, one can expect that for / = li, the vertex (0,0) remains 
irregular, but now it is not singular, i.e., w(0,0~) 7^ (under a natural non-orthogonality 
condition) is finite. Same happens for other roots A& of Ao(Z) (provided that the rest of 
eigenvalues are "ordered", which has not been proved for arbitrary I > 0). 



6.4. Boundary layer and branching at I = +00. The phenomenon of the boundary 
layer (BL) occurring as I — > 00 is explained by Figured It follows that already at / = 20, 
the first eigenfunction ^(y) of B* has a matching structure of the Boundary Layer for 
y m I with the constant eigenfunction ipoiu) = 1 ( see (15.191) ). which dominates for / < 12. 

The BL structure here is simple and standard; see more details in Section 17.21 where 
the BL is constructed for the parabolic PDE. Briefly, asymptotically sharp, in the BL for 
y ~ I 3> 1, equation (16.41) for |A o (0l <S 1 being asymptotically small (see (16.191) ) reads 

-(**)(*) -I/(**)' + ... = W = V(z) + ..., z = l*y, where 

(6.13) 

-vW-\v = o, ze(o,Z3) ; v = v = 0, z = n- 7(0) = 1. 
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First eigenfunction of B : 1=20, X Q (20)=0. 0000041 47.. 



1.2 - 



¥(y) 



0.8 



0.6 



0.4 - 



0.2 



Inner Region 



Boundary Layer 



v„(y)-v (y) = i 



Matching Region 




5 10 

y 

FIGURE 3. The boundary layer structure of ^(y) for I = 20. 
The BL-function V is given explicitly (cf. (I7.10p ): denoting b = 2~ 5 / 3 and a = V3b, 

sin(a/ 4 / 3 )] 



,i 4 / 3 r 



V(z) = C 3 - e" fez [d cos{az) + C 2 sin(az)] , C 
(6.14) 4/3 

e fr' 4/3 [ s in(a; 4 / 3 )+-)=co S (aZ 4 / 3 )] 



<■'" |cos(a/ 4 / 3 ) ^= siu(. 



1-e 



M 4 / 3 



^3 

cos^ 4 / 3 )--^ sin(aZ 4 / 3 )] : 

1 



1-e' 



1-e' 



6(4/3 [cos(a/ 4 / 3 ) L sin(a/ 4 / 3 )]' 



(6.15) 



[cos(ai 4 / 3 )— ^= sin(ai 4 / 3 )] j.— c |u»^w ■ ; — 

The eigenvalue Ao(/) for I ^> 1 is then obtained by matching in R: extending = 
^qH(1 — y), where H is the Heaviside function, we have (see Section E3] for details) 

(n) (4) (y) = (n) {4) (y)H - (nni)5(y - o - wwfo-o, 

where (#5)"(0 = ^"(Js) + = I V'"0) + ... . 

Then (JS3D for A (/) reads 

(6.16) B*^*-/V ,/, (/t)5(y-/)-/iv"(/s)5 , (z/-0 + - = A (0^ in M - 
Finally, substituting = ipQ + w = 1 + w, where W-Lipo, yields 

(6.17) B*w - I V'"(l?)5(y -I)- l*V"(l*)6'(y - I) + ... = X (l) + ... . 

Hence, by the orthogonality condition to i^o(y) = F(y) (recall that B*m;_L?/>o), (16.171) 
implies the following asymptotic expression for the eigenvalue: 

Ao(i) = -{l V"'(ll)5(y -I)- liV"0)8'(y - I), A) + - 



(6.18) 



-i v'"{u)^{i) + iw\uw Q {i) + 
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Overall, using the asymptotics (15.121) . (15 . 1 3[) for the rescaled kernel and taking into account 
(16.141) . we obtain from (16.181) the following typical approximate oscillatory behaviour: 

(6.19) A o (0 ~ Jse"* 14/3 cos (So /s) for / > 1, d ~ d + b, b ~ ^±2, 

where we omit constants. In fact, (16 . 191) shows how an /-branch of the first eigenvalue 
Xq(1) bifurcates at I = +oo from A (oo) = for the operator B* with the spectrum (15. 15)) . 
m = 2. It is not that difficult to show that similar branching occurs from any eigenvalue 

(6.20) A fc = -f for k = 0,1,2,... =s> 3 X k (l) for I > 1, X k (oo) = -f, 

where ( I6.19P is evidently the only eigenvalue that can change sign. The branching happens 
in the framework of classic perturbation theory for linear operators (see Kato [39]), so we 
do not treat this any further. Actually, the branching (16.20)) or (16.101) is the only for us 
reason for the operator ( 16. If) to admit many real eigenvalues for various values o£j I > 0. 

7. B I- HARMONIC PDE FOR m = 2: BLOW-UP ASYMPTOTIC DERIVATION OF 

Petrovskii-type regularity criterion by EIGENFUNCTION EXPANSION 

Regardless nonexistence of a definite answer concerning regularity of the vertex of 
the fundamental parabolae (14.131) . nevertheless, we next show that there is a way to 
move towards a Petrovskii-like criterion for the case of more expanding lateral boundaries 
for functions (14.71) . However, the results in Section [6] (treated according to the limit 
I — 1(t) — > +oo) imply, in view of the oscillatory behaviour of the first eigenvalue in 
(16.191) . that a direct regularity conclusion for a given <p(r) is not possible in principle. To 
get the regular vertex (0,0), a procedure of oscillatory cut off of the boundary must be 
performed. 

At this moment, in Appendix A, we continue our study with a certain brief, rather 
questionable, and even controversial suggestion, which nevertheless is quite attractive 
and actually goes along the lines of Petrovskii's barrier ideas, so we fairly believe this is 
a right (and unique) place for such a discussion. 

As a next step, we return to the spectral methods that cover both regular /irregular 
issues for ( 14. ip . which indeed are more complicated than for the heat equation (I3.ip . 

7.1. Two-region expansion. Thus, we return to the rescaled problem (14.121) for m = 2. 
As usual in any matching asymptotic analysis, this blow-up problem is solved by matching 
of expansions in two regions: 

(i) Inner Region, which is situated around the origin y = 0, and 

(ii) Boundary Region close to the boundaries y = ±<^(r), where a boundary layer occurs. 

Actually, such a two-region structure, with the asymptotics specified below, defines the 
class of generic solutions under consideration according to (14.31) (those in (14.41) are all 
different). We begin with the simpler analysis in the Boundary Region (ii). 

16 We then state an open problem: is there any direct proof that (|6.1|) has real eigenvalues only, and 
the eigenfunctions form a complete and closed set in L 2 (7j)? Or this is hopeless and branching together 
with global extensions of the branches (also an open problem in general) is the only reason. 
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7.2. Boundary Layer (BL) structure. Sufficiently close to the lateral boundary of Qq, 
it is natural to introduces the variables 

(7.1) z = 7$t) and r ) = w ( z i T ) =^ ^ = -^r - Jzw* + ^ 

We next introduce the BL-variables 

(7-2) f = <pt(r)(l -z), y^(Y) d T- = ds, an d r) = p(s)g(£, s), 

where p(s) is an unknown slow decaying (in the same natural sense, associated with (14.81) ) 
time-factor depending on the function ip(r), e.g., ~ as a clue. On substitution into 

the PDE in (17. ip . we obtain the following perturbed equation: 

9s = Ag - \ <p-*£gs - f (l - &r*)# 

V / 1 ' 

-| V'^vUm ~^9, where Ag = -g^ + \ g' . 

As usual in boundary layer theorjQ, following (14.31) . we are looking for a generic pattern 
of the behaviour described by (17. 3p on compact subsets near the lateral boundary, 

(7.4) |£| = o(</?~5( T )) \ z - 1| = o(<^~t(r)) as r — > +oo. 

On these space-time compact subsets, the second term on the right-hand side of (17.31) 
becomes asymptotically small, while all the others are much smaller in view of the slow 
growth/decay assumptions such as (14. 8p for ^(r) and p(s). 
Then posing the asymptotic behaviour at infinity: 

(7.5) g(£, s) — > 1 as £ — ► +00 (Hypothesis I for generic patterns), 

where all the derivatives are assumed to vanish, we arrive at the problem of passing to the 
limit as s — > +00 in the problem (17. 3p . (17. 5p . Assuming that, by the definition in (17. 2p . 
the rescaled orbit {g(s), s > 0} is uniformly bounded, by classic parabolic theory [19] , 
one can pass to the limit in (17. 3p along a sequence {s^} — > +00. Namely, by the above, 
we have that, uniformly on compact subsets defined in (17.41) . as k — > 00, 

(7.6) g(sk + s) — > /z(s), where /i s = A/i, h = = at £ = 0, /i|^ =+OC) = 1. 
The Zimit (at s = +00) equation obtained from (17. 3p : 

(7.7) h s = Ah = -h^K + \h(. 

is a standard linear parabolic PDE in the unbounded domain K. + , though it is governed 
by a non self-adjoint operator A. We need to show that, in an appropriate weighted 



It was always key for PDE theory. Shortly after the Blasius construction (1908) of the exact self- 
similar solution for the two-dimensional boundary layer equations proposed by Prandtl in 1904 (see 
references in [32| p. 48]), similarity solutions of linear and nonlinear boundary- value problems became 
more and more common in the literature. 
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Boundary Layer, m=2: the first stationary pattern g {%) 




FIGURE 4. The unique stationary solution go(£) of the problem (|7.8p . 

L 2 -space if necessary, under the hypothesis (17.51) . the stabilization holds, i.e., the w-limit 
set of {h(s)} consists of equilibria: as s — ► +00, 

(h(t s) -> g Q (£), where Ag = for £ > 0, 

(7.8) < 

[00 = 00 = for £ = 0, g (+oo) = 1. 

The characteristic equation for the linear operator A yields 

(7.9) -A 4 + ±A = Ai = and A 2)3 = ^ ( - | ± i 
This gives the unique solution of (I7.8p . shown in Figure HI 

(7.10) g (0 = 1 - e~^ [cos (0) + j= sin (0)} . 

We do not concentrate on this stabilization problem ( 17.81) . which reduces to a standard 
spectral study of A in a weighted space. Let us mention that, in the metric of L 2 (M + ), 
the operator has a clear exponential stability feature. Consider the eigenvalue problem 

(7.11) Aip = -^ (4) + \tf = Xip- 

As customary, multiplying this by ip, next taking the complex conjugate of the equation 
and multiplying by ip an d summing up the resulting equalities yields 

(7.12) (a + a) / h 2 = -2j m 2 + \ jw + m = -2 / \r\ 2 < 0, 

i.e., Re A < for any A G er(A) (this does not take into account a continuous spectrum). 

Actually, the convergence (17.61) and (17. 8p for the perturbed dynamical system (17.31) 
is the main Hypothesis (H), which characterizes the class of generic patterns under 
consideration, and then (17. 5p is its partial consequence. Note that the uniform stability 
of the stationary point go in the limit autonomous system (17.71) in a suitable metric will 
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guarantee that the asymptotically small perturbations do not affect the omega-limit set; 
see [Ml Ch. 1]. 

We stop further discussion concerning the passage to the limit in ( 17. 3j) and summarize 
the conclusions as follows: 

Proposition 7.1. Under the given hypothesis and conditions, the problem ( \7.3\\ admits a 
family of solutions (called generic) that satisfy (17.81) . 

We must admit that such a definition of generic patterns looks rather non-constructive, 
which is unavoidable for higher-order parabolic PDEs without positivity issues. Note 
that, according to (I3.14p . a similar obscure issue appears even for the heat equation in 
(13. ip . One can expect that (17. 8p occurs for "almost all" solutions, excluding just those 
that have eventually the faster vanishing first Fourier coefficient that some others. Then 
another BL is needed, but we do not intent to describe an invariant manifold structure of 
such a thin solution set (recall that (17. 3p is a difficult non- autonomous PDE). 

7.3. Inner Region analysis: towards the regularity criterion. In Inner Region, we 
deal with the original rescaled problem (I4.12p . Without loss of generality, for simplicity 
of key calculations, we consider symmetric solutions defined for y > by assuming the 
symmetry conditions: 

(7.13) Vy = Vyyy =0 at ?/ = 0. 

In order to apply the standard eigenfunction expansion techniques by using the or- 
thonormal set of polynomial eigenf unctions of B* given in (I5.19p . as customary in classic 
PDE theory, we extend v(y, r) by for y > (f(r): 

,7,4) Hv ,r) =, to ,r W »Kr) -») = {"(». J) £ ° y f<^> 

where H is the Heaviside function. Since v = v y = on the lateral boundary {y = (p(r)}, 
one can check that, in the sense of distributions, 

V T = V T H, Vy = VyH, Vyy = VyyH, 

(7.15) ^yy = Vyyy H -Vyy\ y= ^5{y -If), 

%vm = v yyyy H ~ v yyy \ y= ^{y - <p) - v yy \ y= J'(y - <p). 
Therefore, v satisfies the following equation: 

(7.16) v T = B*v-Vyyy\ y=ip 5{y-v)-Vyy\ y=v 5\y-y) in l + xl + . 

Since, obviously, the extended solution orbit (17.141) is uniformly bounded in L 2 p * (R), we can 
use the converging in the mean (and uniformly on compact subsets in y) the eigenfunction 
expansion via the generalized Hermite polynomials (15.191) : 

(7.17) v(y,r)= £ a k (r)r k (y). 

(fc>0) 
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Substituting (17.171) into (17.161) and using the orthonormality property (I5.20p yields the 
following dynamical system for the expansion coefficients: for all k = 0, 1, 2, ... , 

(7.18) a' k = \ k a k ~ Vyyy\ y=tpiT) (5(y - <p(T)),ip k ) - v yy \ y=tp(T} (6'(y - (p),ip k ), 

where = — | are real eigenvalues (15.151) . Recall that < for all k > 1. More 
importantly, the corresponding eigenfunctions ip k {y) are unbounded and not monotone 
for k > 1 according to (I5.2ip . Therefore, regardless proper asymptotics given by (I7.18p . 
these inner patterns cannot be matched with the BL-behaviour such as (17.51) . and demand 
other matching theory (since these are not generic, the latter is not developed). 

Thus, bearing again in mind (14.31) . one needs to concentrate on the "maximal" first 
Fourier generic pattern associated with 

(7.19) k = 0: A = and ^ (y) = 1 {Mv) = F(v))- 

Actually, this corresponds to the centre subspace behaviour for the equation (17.181) : 

(7.20) v(y,r) = ao(r) +w ± (y,T), where w ± G Span{^, k > 1}, 

and w _L (y,r) is then negligible relative to o>o{t). This is another characterization of our 
class of generic patterns, Hypothesis II. The equation for ao( r ) then takes the form: 

( 7 - 21 ) a 'o = - v yyy\ y=ip{r) M^)) + v yy\ y= ^ T) %^i. T ))- 

We now return to BL theory established the boundary behaviour (I7.2p for r ^> 1, which 
for convenience we state again: in the rescaled sense, on the given compact subsets, 

(7-22) v(y, r) = p(s)g ( V l(r)(l - ^)) + ... . 

By the matching of both Regions, one concludes that, for such generic patterns, 

(7.23) ^^1 as r^oo =^ p(s) = a (r)(l + o(l)). 

Then the convergence (I7.22p . which by regularity is also true for the spatial derivatives, 
yields, in the natural rescaled sense, 

»wL w P( s )¥^(7)7i -> a (r)<fl(r)^i, 71 = ^'(0) > 0, 

(7.24) y ~ v( ' 

Vyyy\ y=tp{T) -> -p(s)(f(r)j 2 ~> -flo^M^, 7 2 = ^o'( )- 

Eventually, this leads to the following asymptotic ODE for the first expansion coefficient 
for generic pattern^ 

(7.25) J = G 2 (^(t)) = 7 2^(t)V'o(^(t))+7i^Wo(^)) + --- for r » 1. 

For instance, this gives a first easy condition of regularity of (0,0): for the generic 
patterns, this is the "negative" divergence of the following integral: 

00 

(7.26) f G 2 ((p(r)) dr diverges to —00 =^> a (r) — > as r — *• +00. 



18 This result can be stated as a theorem: For the prescribed above class of generic solutions, the first 
Fourier coefficient satisfies... ; as we have mentioned, we avoid such non-constructive, but rigorous, ones. 
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In order to obtain more practical and sufficiently sharp conditions of regularity, we will 
use the expansion (15.121) of the first eigenfunction ipo(y) = F(y), which on substitution 
into the right-hand side in (17. 25ft . where both terms are equivalent, yields 

(7.27) j£ =^{r)C 3 cos(b oV l{T) + C 4 )e- d ^ 4/3 ^ + ... for r > 1, 

with some 7 ^ and constants depending in an obvious way on C\ t 2 in (15 . 1 2[) and 
other parameters from (15.131) . Integrating implies that 

(7.28) \n\a {T)\=^j^{s)C lC os(bo^{s) + C 2 )e' d ^ 4/3 ^ds + ... for r > 1. 

Irregularity condition. This is straightforward: (I7.28P implies that the limit of ao(r) as 
t — > +00 can be arbitrary (i.e., not necessarily zero), if the integral converges. Therefore, 
up to the achieved accuracy of the expansions and matching, we state the following 
condition of the irregularity of (0,0) (all the irrelevant constants are omitted): 



(7.29) 



2 4 

(0,0) is irregular if J <p~3 (s) cos (60V 91 e~ dov ( s 'ds converges. 



Then, obviously, w(0,0~) 7^ (a natural non-orthogonality condition required) is finite, 
i.e., the vertex (0,0) is irregular, but non-singular. This still looks pretty similar to 
Petrovskii's criterion in (13. lip . 

Regularity condition: oscillatory cut-off of ^{t). Thus, assume that f{r) is close to the 
desired (and still unknown) critical irregularity /regularity situation, and that we follow 
the first Fourier coefficient equation (I7.27P or (17.281) . within the given accuracy. Then, 
for using (17.281) for the purpose of the regularity analysis, we observe that the presence 
in the integral the oscillatory factor can always cos(-) violate the (uniform) divergence 
to —00 hypothesis in (I7.26p . Indeed, then In |ao(r)| — > ±00 along different subsequences 
{t = Tj^} — > +00, and the limit to +00 does not guarantee (I7.26P (of course, for such 
limits, our matching with the BL can be violated as well, but we do not discuss this issue 
here). In other words, in the case of oscillatory kernels, a pure divergence of the integral 
in (I7.28P does not guarantee the point regularity. This is the principal difference with the 
positive kernel case; cf. (I7.36P below with no oscillatory component in the integral. Let 
us note that, in general, the standard divergence of (I7.28P will ensure that that the vertex 
(0, 0~) is irregular and is also singular in the sense of (14. 15[) . Moreover, u(x, t) has infinite 
oscillatory behaviour as (x,t) — > (0,0~), i.e., eventually oscillations get both ±00. 

Thus, under the accepted hypothesis, the oscillatory part of the rescaled fundamental 
kernel (I5.12p of the operator under consideration can violate the regularity of boundary 
point for any slow growing factor <£>(t) satisfying (I4.7l) - (l4.9p . even if the integral diverges. 
As we have seen in Section [6j for higher-order operators, the regularity condition "above" 
the fundamental parabolae (14.131) becomes rather subtle and sensitive. Indeed, accord- 
ing to the integral in (I7.28p . arbitrarily small perturbations of v?(t) can "switch over" 
regularity to the irregularity and vice versa. 

To guarantee regularity, an extra "oscillatory cut-off" of <p{r) is then necessary to be 
introduced after a necessary remark. 
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7.4. Remark 1: an analogy with elliptic theory. This principal difficulty concerning 
the kernels of changing sign has a known counterpart in regularity elliptic theory for (12. ip . 
Namely, it was first shown in [70] (1986) for m = 2, N > 8 that the vertex of a cone can 
be irregular (singular) if the fundamental solution of L(d) for N > 2?d0 

(7.30) F{x) = ^(^)\x\ 2m - N , xeR N \{0} (N>2m) 

changes sign; see further comments and references in [B"8"| § 1]. Then in elliptic theory, the 
regularity analysis is performed for a restricted subclass of operators L(d) called positive 
with weight F. This actually means that F > [SSJ § 3] (cf. a positivity-like condition 
in [27] for m = 2 and [21J), so, until now, refined elliptic regularity results for kernels 
of essentially (e.g., infinitely many) changing sign seem be unavailable. As in classic 
theory [491 [71], this analysis demands constructing corresponding generalized Hermite 
polynomials as eigenfunctions of an adjoint pencil of linear operators along the lines of 
those obtained for hyperbolic equations; see (32[ p. 254] and [29] (for u tt = —u xxxx \ see 
the end of Section [S] below). The oscillatory cut-off of the boundary will then be essential. 
An example of such an elliptic evolution approach for studying singularities of (12.21) in 
K 2 , q G (0, 1) (with many additional references presented) is given in |10j . 

Of course, it is well-known and this is a classic matter, that general theory of operator 
pencils associated with corner singularities has been well-developed for a number of elliptic 
equations; see [521 [53] and [27J for references to other papers and related monographs and 
as a source of further extensions to be traced out by the MathSciNet. For other types 
of parabolic or elliptic problems mentioned above, new operator pencils appear. In the 
parabolic case, we do not intend to change the operator — D\ in related lineo, but then 
need instead to perform a refined oscillatory cut-off of the lateral boundary. 

7.5. Optimal regularity conditions for generic patterns. Thus, to guarantee the 
regularity, we perform the procedure of oscillatory cut-off of a given function f(r). 
Namely, this means that a smooth ^(r) replaces fir), for which 

(7.31) 7 C a cos (6o0(r) + C 2 ) « [*fC x cos (b <pl(r) + C 2 )]_ 

asymptotically sharp up to absolutely convergent perturbations in the integral (I7.28p . 
where [•]_ denotes the negative part. This is necessary to cut-off the positive part of 
the diverging integral. Figure [5] schematically explains how (17.311) works (for 7C1 > 0) 
by cutting off all positive waves of the cos-function for r 3> 1 and creating "almost 
discontinuous" (jumping-like) <p(r), which in the figure is denoted by <p cu t( T )- Using 
necessary smoothing at the points, where cos(-) changes sign, for any such monotone 
increasing ip(r), the corresponding <£>(t) can be also chosen increasing. In other words, the 

19 For N — 2m, the structure of the fundamental solution is different and more signed-determined: 

F(x) = «ln^ +*(]§t) (k(N) = const.), 

that makes the case special admitting arbitrary operators L(d) 68, § 8]; Wiener's m = 1, N = 2 included. 

20 Obviously, unlike elliptic theory, this is not possible: no higher-order parabolic equations with smooth 
coefficients can have positive kernels (otherwise, comparison and the MP would be inherited). 
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Constructin oscillatory cut-off of $(t) 




FIGURE 5. A schematic view of a proper oscillatory cut-off of <p(r). 

resulting <f(r) sufficiently fast jumps over those intervals of the length /, which in Section 
[6] were determined as leading to the irregular vertex. We do not pay special attention 
to the question on how such smoothly jumping boundary can affect the boundary layer 
behaviour described by the equation (17.31) . where the term f'/f may then play a role 
(fortunately, not dominant). However, it is clear that, since the correct intervals of the 
behaviour of (p(r) get arbitrarily long as r — ► oo, we have enough time for establishing 
the BL-structure, so that formulae (18.51) remain true closer to the end of each of them, 
and further analysis applies. In addition, it is also clear that, once the solution gets very 
small by (I7.25P at the end of a fixed such interval in r, a smooth monotone increasing 
jump of (p(r) at the end point cannot essentially affect the value of the solution during 
such a short interval of time. This means that, on the next good interval, the solution 
remains small and continues to be governed by the eigenspace Span {ip^}- 

Thus, according to the oscillatory cut-off (17.311) . the condition of regularity is as follows: 



(7.32) 



(0, 0) is regular if J (ps (s) cos (&0V 93 ( s ) + C2) e do<p ^ ds diverges 



There is a clear gap between the irregularity (I7.29f) and the regularity ( 17. 32ft conditions, 
so that an Osgood-Dini-like criterion cannot be available in principle. 

With such a procedure, similar to Petrovskii's presentation ( 13.91) of the boundary reg- 
ularity results, we deduce that: 



(7.33) 



(i) R(t) = 


3-1 2^ ( 


-*)* [In | ln(— ^ = 


=>■ (0, 0) is rej 


jular, and 


(ii) R(t) = 


(3-1 2^ 


+ e) (-*)* [ln|ln(-t) 


}~\ e>0 


(0,0) is irregular, 
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where, in the regular R(t), the oscillatory cut-off meaning (17.311) must be assumed. 



7.6. Remark 2: on other asymptotic patterns and related regularity. It follows 
from the dynamical system (17.181) . that, in general, there exist other (non-generic) as- 
ymptotic patterns corresponding to the behaviour on each of a ID stable eigenspace of 
B*, which hence are not governed by the expansion (I7.20p . Such a behaviour then will 
generate its own regularity/irregularity conditions with or without cut-offs of the cor- 
responding similar integrals. Indeed, those asymptotic patterns will demand a different 
"BL-like" theory, for which (I7.5p is not true. Since these are not generic and belong to 
the subspace of co-dimension one, a correct posing of the IBVP in Q for such regular 
solutions is quite tricky and demands a priori unknown and non-constructive conditions 
on initial data uq, so we have no reasons to take those into account. On the other hand, 
describing countable sets of various blow-up singularities is a serious problem of modern 
PDE theory, but currently this has a little to do with the boundary regularity analysis. 

7.7. Remark 3: for the heat equation, the spectral— BL blow-up approach is 
sharp. Let us very briefly repeat and list the main steps of the above blow-up analysis 
for the classic problem (13. ip to prove (14.51) . Thus, (14 .6p reads 

R(t) = H)^(r), 
where, in (I4.12p . we get the classic Hermite operator [TT1 p. 48] 

W = Dl-\yD y in Lj.(R), p*(y) = e - y 4, a (B*) = {\ k = -f, k = 0, 1, 2, ...}, 
with X>(B*) = Hp,(M), etc. Similar to (17.21) . the BL variables and asymptotics are now 
£ = <p 2 (l-z) = <p(<p-y) and w(z, r) = ao(r)g(£, r) + ... , 

(7.34) 

where g(£, r) -> g (£), r -> oo and g (£) = 1 - e ^ 2 [g$ + \g' Q = 0). 
Using two formulae in the first line in (17.151) with 

Vyy = VyyH -Vy\ y = ^5{y -<f), 

we arrive at the equation (the analogy of (17.161) ) 

v T = B*v + v y \ y=ip S(y-(p). 

Hence, asymptotically via the BL structure (17.341) . for the generic patterns as usual, the 
analogy of (18.61) takes much simpler form: 

(7-35) J = - 7l ^oM + ... =► % = -&V*-£ + ~, 

where 71 = ^(0) = | and ipo(y) = F(y) = 2^ e ~ y2 ^ ls ^ e positive rescaled Gaussian 
kernel of the heat operator. Thus, (I7.35P justifies that, for the generic patterns, the 
regularity criterion reads 



(7.36) 



(0, 0) is regular iff J ip(r) e v 4 dr = +00 
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(no cut-off is necessary in this non-oscillatory case). One can see that (I7.36P is equivalent 
to that in (13. lip or (I3.12p . Recall that there exist other asymptotic patterns related to 
the stable subspace of B* of co-dimension 1, which, as usual, is not taken into account. 



8. ON EXTENSIONS TO 2mTH-ORDER POLY-HARMONIC AND OTHER OPERATORS 

We discuss the extensions of the asymptotic method to other PDEs. 

8.1. Poly-harmonic equation. For (11.41) . the spectral properties of the operator 

(8.1) B' = -(-A)»-l S .V 

are given in [TT] . A sharp asymptotic expansion (similar to ( 15.121) ) of the radially sym- 
metric rescaled kernel satisfying 

(8.2) BF = in R N , JF=1 Mv) = dPf (v), |/?| = 0, 1, 2, ... , 

where (3 is a multiindex in R N , is also not that difficult to get. Then ( 18.21) implies that 
the adjoint basis {ipt} consists of generalized Hermite polynomials. 
In the radial case, the critical R(t) for boundary regularity is about 

1 2m— 1 _ JL, 

(8.3) R(t) = 0* (-t)^ [ In | ln(-t)| ] 2m , where C* = d a > is a constant, 

and the oscillatory cut-off is assumed for the regularity. As usual, replacing C* by C* + e 
for any e > makes (0,0) irregular. 

Some calculations become more involved, especially in the case of non-radial Qo, where 
the spatial shape of the shrinking cusp will affect the right-hand side of (17.271) . where the 
BL structure, though having similar variables such as (17.21) . is now governed by compli- 
cated elliptic problems instead of the ODE ( 17. lip . For instance, if the boundary of Qo is 
given by the characteristic fundamental paraboloid in the x and y variables, respectively 
(i.e., tp(r) = 1): 

( 8 - 4 ) £(M=2m) a « X ° = (-*) E(\ a \=2m) a »y a = 1- 

The case a a > for any a corresponds to Mihailov's [721 ESI El] and Kondrat'ev's [IB] 
cases. Hence, the regularity in the strong sense of the vertex (0, 0) depends on the 
spectrum of the operator (18. ip in the domain I a with the boundary given by (18.41) . Similar 
to the analysis in Section [BJ we have that, by Poincare's inequality, (0, 0) is regular if 
the diameter of I a is not that large (since the eigenvalues of B* are close to those for 
— (— A) m < 0). On the contrary, since B* is not self-adjoint and sign-definite, there are 
coefficients {a a > 0}, for which the vertex (0,0) is not regular in the classic sense. 

If some of the coefficients {a a } are negative (this is Feigin's case |25j), B* is then posed 
in an unbounded domain with a non-compact boundary, and its spectral theory becomes 
more involved. Nevertheless, for such sufficiently "thin" domains, we may expect the 
vertex be regular, and sometimes irregular otherwise. 
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8.2. Linear dispersion equations. There is no much difference in studying the bound- 
ary regularity for odd-order linear PDEs. For instance, as a simple example, consider the 
third-order linear dispersion equation in a similar setting: 

(8.5) u t = u xxx in Q = {-R_(t) < x < R + (t), -1 < t < 0}, 

where R±(t) are positive smooth functions on [—1, 0), continuous on [—1, 0], and R±(0) = 
0. Since the operator is anisotropic, the functions R±(t) are essentially different. 
According to odd-order PDE theory (see e.g., [21]), for (18.51) . it is allowed to put the 
following Dirichlet boundary conditions: 

(8.6) it = u x = at x = R-(t) and u = at x = R+(t), i G [-1,0). 

Since, as we will show, the main singular phenomena occur at the right-hand "oscillatory" 
lateral boundary {x = R + (t)}, we concentrate on this analysis, and, in general, can put 

(8.7) R_(t) = -oo, 

so excluding this part of the boundary from consideration. Then, the regularity criterion 
of (0, 0) will solely depend on the behaviour of R+(t) as t — > CT. 

Thus, taking into account the right-hand lateral boundary, as in (14.61) . we perform the 
first rescaling: 

to ox R +® = ( _t )^+( r )' y = rfvs ' r = - ln (- £ )> M ( x > *) = v (y> r )> 

where v T = B*t> = v yyy — | yv y , v = at y — y+(r), r > 0. 

Hence, unlike parabolic equations, for (18.51) . another (but indeed similar in many places) 
version of Hermitian spectral theory for B* is necessary to get optimal behaviour of R+(t) 
as t — » _ for the regularity/irregularity of (0, 0). 

First, as above, the regularity criterion is associated with the structure of the funda- 
mental solution of (18. 5p that has the similarity form 

(8.9) b(x,t)=t-$F(y), y=tfe, 
where F = Ai(y) is the Airy function satisfying 

(8.10) BF = F'" + | (yF)' = in E, J F = 1 F" + | yF = 0. 

Hence, the rescaled kernel F(y) has an exponential decay as y — ► — oo only, and as 
y — ► +oo is oscillatory according to the asymptotics 

dlyl-^e-^ 3 ' 2 + ... as y -> -oo, 

jT3 [C 2 cos(do2/^) + C3 sin(do^)] + ••• as y -> +00, 

where c? = and are some constants of non-essential values. 

Secondly, this (see extra details in [28], § 9.2] and [SB] ) generates very oscillatory as 
y — ► +00 and unbounded (for any k > 1) eigenfunctions of B in (18.101) (cf. (I5.16P ) 

(8.12) ^ k ( y ) = t^D k y F(y) (F — M(y)), with a(B) = { - |, k = 0, 1, 2, ...}. 
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(8.11) F(y) 



For existence of such a discrete spectrum, a proper "radiation" condition at y — +00 is 
posed, which excludes "non-oscillatory" (polynomial-rational) bundles that are different 
from those in ( 18. lip . The "adjoint" linear operator 

(8.13) B* = D 3 y - \ yD y , with *{B*) = a(B) = { - §, fc = 0, 1, 2, ...}, 

has a complete set of eigenfunctions {ipt} that are generalized Hermite polynomials. In 
general, then the bi-orthonormality property (15.201) demands using Hahn-Banach Theo- 
rem on extensions of continuous linear functionals (or similar techniques of regularization 
of oscillatory integrals in a v.p. or c.r. sense) and other constructions (this is not that 
crucial for in what follows). Both operators are defined in special weighted L 2 spaces. 
One can see that B* is not adjoint to B in the standard metric of L 2 , since, obviously, 
then B* = — B + ~ /, so B — 1 1 is skew-symmetric. In fact, B* is adjoint to B in a space 
with the indefinite metric given by 

(8.14) (v, w)* = fv(y)w(-y) dy. 

This case corresponds to the decomposable space with indefinite metric with straightfor- 
ward majorizing one and is treated as rather trivial; see Azizov-Iokhvidov |7J for linear 
operators theory in spaces with indefinite metric (the theory was initiated in the 1940s 
and 50s by Pontryagin and Krein). Then the domain of B* is defined as H%, etc. 

Thus, it is key that the bounded operator B* : H 3 , — > L 2 , admits a complete set of 
polynomial eigenfunctions $* = {ipl {y)}, which are constructed similarly to those obtained 
in ( 15.191) . Quite analogously, these polynomials are used in our regularity analysis. 

Boundary Layer at y = R + (t). Similar to ( 17. ip . ( 17. 2p . the BL- variables are 

(8-15) z = ^y e = 4W(l-4 w(z,T)=p(r)g^,r). 

This yields the asymptotic problem of stabilization to a unique profile go (£) (as usual, we 
omit this not that easy analysis), where, as seen from (I8.8p . 

(8.16) tf» - § g' = 0, «7o(0) = 1, g (+oc) = 1 g {£) = 1 - e"^. 

Inner Region: the regularity criterion for R+(t). Using the formulae as in ( 17. 15ft . we ob- 
tain here the equation 

(8.17) v T = B*v + v yy \ y=v+ 5(y-<p + )+v y \ y=v+ 5'(y-<p + ) in R x R + , 

and use the eigenfunction expansion ( I7.17P via the corresponding generalized Hermite 
polynomials that leads to the dynamical system such as ( 17.1 8ft . Looking for the generic 
patterns with the representation (I7.20p . (18.171) yields 

(8.18) a' = v yy \ y=ip+{T) Mv+( r )) - ^I^ +(t) ^o(^+( t )) (Mv) = F (v))- 

By BL theory and matching, we can use the fact that, in the rescaled sense, on the given 
compact subsets, with first derivatives 

(8.19) v(y, t) = a (r)g (/Ur)(l - ^)) + ... . 
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Eventually, this leads to the following asymptotic ODE for the first expansion coefficient 
for generic patterns: 

(8.20) J = G 3 (cp + (r)) = 7 2^ + (r)^o(v + (r)) + 7i 4W^(v+(t)) + •■■ for r » 1. 
Since by the asymptotics (18. lip both terms are similar, this yields 

(8.21) & = 7<pl{T)C4COB(do<pl(T) + C B ) + ... for r > 1. 
This leads to the following condition of irregularity of (0, 0): 



(8.22) 



— — 

(0,0) is irregular if J y?+(s) cos (d (p+(s)) ds converges. 



There is no log-log factor in the critical R+(t) and just a single log occurs: in view of the 
influence of the oscillatory factor cos(-), the point (0,0) 



8.23) for R+(t) = (— t)s | ln(— 1)| 7 : regular for 7 < | and is irregular for 7 > 



where, for the regular case, the corresponding oscillatory cut-off on <£>(t) according to the 
rule as in (17.3 ip is assumed, i.e., then R+(t) reads R+(t). 

Concerning the left-hand lateral boundary {x = —R-(t), t E [—1,0)} (setting R+(t) 
■ x : . the standard exponential decay behaviour of the kernel (18.111) as y — > —00 is 
supposed to give a result which is similar to the heat equation Gaussian, with F(y) ~ 
e -?/2 / 4 . Since the decay rate is slower in (18. lip governed by the power \y\^ , bearing in mind 
the integral criterion in (I7.29p . one can derive the following conclusion on the regularity: 



(8.24) f°° ... e - do ^-^ 3/2 ds = 00 it_(t) = (^) ~ 3 (-t) I [in | ln(-t) 

where the constant (^^) 2/ ' 3 cannot be increased without losing the regularity of (0,0). 
Since the kernel F(y) is of constant sign in this limit, no oscillatory cut-off is necessary. 
Thus, the third-order PDE (18. 5p is the last one, for which for a (one-sided) regularity 
criterion, there occurs a non-oscillatory kernel and the result looks similar to the heat 
equation one (13. 9p (though no sub-solutions like (I3.5P can be used, and, hypothetically, a 
majorizing order-preserving flow as in Appendix A can be discussed). As we have seen, 
for the bi-harmonic equation (14.11) . both lateral boundaries, in general, are assumed to 
undergo oscillatory cut-offs for the regularity. The same is true for the fifth-order linear 
dispersion equation and the tri-harmonic one 

(8.25) Uf u xxxxx and v,f u xxxxxx , etc., 

on their and others oscillatory properties, see [321 Ch. 4 and 3]. 

More complicated spectral theory of multi-dimensional odd-order operators B* and B 
is necessary for the study of the boundary regularity phenomena for the iV-dimensional 
counterparts of linear dispersion PDEs such as 

(8.26) u t = (Au) xi + ... or u t = -(A 2 u) Xl + ... in Q . 
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8.3. On a quasilinear fourth-order diffusion equation. Here we briefly explain nec- 
essary changes towards the regularity analysis for quasilinear diffusion-like PDEs. We 
take the ID cubic 4th-order porous medium equation (the PME-4), as a basic model, 

(8.27) u t = -(u 3 ) xxxx in Q Q , 

where Qq is the same as in fl4.ll) and the Dirichlet boundary conditions now read 

(8.28) u = {u 3 ) x = at x = ±R{t), -1 < t < 0. 

We specially have fixed a fully divergent diffusion model (18.271) with the monotone operator 
in the metric of H~ 2 , for which existence and uniqueness is guaranteed by classic theory; 
see Lions [60j Ch. 1,2]. Using the scaling (14. lip yields 

(8.29) v T = B*(v) = -(v 3 ) yyyy -\yv y , 

where B*(t>) is now nonlinear. After the BL-analysis, we return to (18.291) to see how the 
linear operator B* as in (14. 12ft will be recovered. 

Boundary layer. The approach is similar with the BL variables as in (17.21) . where 

(8-30) £=^g(l-z), w(z,r) = p(r)g^r) (z = fa). 

The BL rescaled equation such as (17.31) is easily recovered with the nonlinear operator 

(8.31) g : A(g) = -(g 3 )^ + \ g' = 0, £ > 0; g = (g 3 )' = at f = 0, g(+oo) = 1. 

The proof of existence for the problem ( 18.311) is not easy, but indeed doable by a standard 
shooting approach, though the uniqueness can cause a serious technical problem. To 
avoid an essential declining to classic ODE theory, which is not the subject of the present 
PDE research, we avoid spending a couple, or possibly a few of next pages with various 
elementary and not that easy related ODE calculus. 

Instead, using standard facilities, we again just check these key properties of existence- 
uniqueness of <7o(0 by using the solver bvp4c of the MatLab with the enhanced accuracy 
and tolerances (16. 9p . The sufficient for us and reliable numerical justification is shown in 
Figure El where the semilinear equation is solved: 

(8.32) G = g 3 : -G® + ^ \G\~sG' = 0, G = G' = at £ = 0, G(+oo) = 1. 

The dotted line shows the "linear" go(£) from Figure HJ so that both linear and nonlinear 
profiles are similar, though correspond to different ODE mathematics. 

In general, stabilization in the full model, which is a quasilinear extension of (17. 3p . is 
rather technical and involved, but a class of such generic solutions forming the necessary 
boundary layer is indeed obtainable. Moreover, instead of the corresponding Proposition 
I7.1[ which is difficult to prove in a sufficient generality, we may view this result as a 
definition of the class of generic solutions under consideration. For this class of solutions, 
the following holds on appropriate compact subsets near the boundary (cf. (I7.22p ): 

(8.33) u ( yjT ) = p ( T ) flb (3^(i__jL y )) + ... as r - +oo. 
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Boundary Layer, m=2: the first pattern GJQ 
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FIGURE 6. The unique stationary solution Gq(£) of the problem (|8.32j) . 

Thus, in the nonlinear case, the BL analysis is rather similar and demands only the 
extra scaling factor in the BL- variable £ in ( I8.30ft ; cf. the linear one ( 17.2ft . 

Inner Region. This is also similar but there is a single nonlinear aspect that we have to 
address to. Thus, the only principal difference between (I8.29P and (I4.12p is that B*(u) is 
a nonlinear operator, so we cannot immediately use the eigenfunctions expansion ( 17. 17ft . 

To this end, we recall that the BL-behaviour ( 18.331) actually implies that, in this critical 
case, we look for the so-called "mesa-like" asymptotics of solutions, which are supposed to 
be almost flat and to be independent of y for r ^> 1 on compact subsets bounded away from 
the BL at \y\ ~ <p{t) (where f!8.33j) takes place). Such mesa-like asymptotics correspond 
to truly critical cases and are well-known in react ion- diffusion theory; see examples in 
[33~t p. 59, 184]. According to this, it is correct to use the following approximation in the 
nonlinear term in B*(i>): 



on compact subsets bounded away from thin domains, where the "non-flat" BL-expansion 
(18.331) is in charge. Then 



so that, to get back to the linear operator B* in the equation, one needs to rescale the 
spatial variable as follows: 





for t > 1 






where 




Then w(r], t) solves the equation that in leading terms coincides with ( 17.161) . 
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where +... denotes higher-order terms due to the sharp approximation (18.341) . In compar- 
ison with the linear case ( 17. 16ft . we observe the additional second term of the order 0{f-j 
for t ^> 1, which is clearly negligible in view of (|4.7p — (j4.10p . Moreover, w v = on the 
subspace spanned by i/jq = 1, for which (I8.34j) and (17.51) hold. 

We next introduce w according to the cut-off as in (IT. 14j) at r\ = <f(r) — > +oo given 
in (18.351) . Then performing it first for the quasilinear equation for w and next using the 
linearization as in (18.361) . we arrive at 

w T = B*w + — \ VWr, 

(8.37) P 

Using in the last terms the BL estimates following from (18.331) (similar to (17.241) ). and, as 
usual, applying (17.231) for the class of generic solutions, the eigenfunction expansion (17.171) 
for (I8.37P yields the following asymptotic ODE for the first Fourier coefficient (a > 0): 

(8.38) a' (r) = 72V^M^M)(^) + 7i4(^ § Wo($=) + - • 

It follow that, unlike a simpler equation (17.251) . due to the extra "nonlinear" scaling 
(18.351) . the Fourier coefficient ao(T) in f!8.38j) is now entering the oscillatory rescaled kernel 
ip and its derivative ip' . Since <p(r) — > +oo and ao( r ) — * for the regularity (see (I7.26P ). 
the first term in (18.381) is leading that yields on substitution the asymptotics (I5.12p . (I5.13P : 

(8.39) a' = 7a Me-^/v^) 4/3 C 3 cos [ b o{^)^ + C i] + - ■ 



The asymptotic analysis of this ODE is similar. For instance, the following rate of con- 
vergence is obtained for the unperturbed "Gaussian" parabola: 



(8.40) (p(r) = l: a (r) ~ 3§ 2-^(lnr)-§ -> 0, 

which shows the obvious result that the fundamental kernel parabolic shape for the shrink- 
ing point of Q is regular. 

In general, taking into account the main exponential term in ( 18.391) only, i.e., fixing the 
simplest approximating model 



8.41) a' = -e-^A/^) 473 (a (r) > 0), 



we see that the critical R(t) corresponds to 

(8.42) R*{t) ~ C(-*)i[ln|ln(-*)|]* as *->0", 

where C > is arbitrary. Of course, unlike the linear case in ( I7.33p . the constant C in 
( 18.421) cannot play a role in view of the scaling invariance of the PDE ( 18.271) . 

(8.43) u = Aii, x = \[Ax {A > 0), 

which changes the constant C without changing the regularity of the point. Similar to 
(I3~T8D - (I3T20D . the scaling (gZSD itself implies (IH^2|) by choosing A(t) = a (r) from (IH^0|) . 
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It is not difficult to perform a more detailed study of the original ODE (18.391) . though a 
too thorough analysis seems excessive. As usual, for the regularity, an oscillatory cut-off 
is necessary. Thus, it turns out that for the quasilinear ( 18. 2 7ft and linear ( 14. 11 fourth-order 
parabolic equations, the regularity conditions (including the concepts and the methods) 
of boundary regularity analysis can be obtained in similar lines. 

An analogous study can be done for the non-divergent thin film equation (TFE-4) 

U X xx)xi 

though its existence-uniqueness theory is less developed, so the results will be more formal. 
Of course, the same techniques apply to the TFE-(2,2), with other distributions of the 
inner and outer derivatives in the differential form 

1H U X x)xxi 

and to other PDEs with more general nonlinearities. Extensions to various quasilinear 
counterparts of the 2mth-order parabolic equations (11.4j) are also possible on the basis of 
the non self-adjoint spectral theory in [17J. 

8.4. Fourth-order hyperbolic equation: regularity via Hermitian spectral the- 
ory for a pencil. We recall that operator pencil theory fully occurred already in the 
regularity study for elliptic PDEs H)] . and was later developed and extended in many 
papers; see [5UI ED ESI EZ] for references. Boundary regularity for hyperbolic PDEs 
is definitely less developed; see last pages of Kondrat'ev-Oleinik's survey [50] of 1083 for 
some references. Note that, for fixed standard characteristic paraboloids posed for (I8.44p . 
the regularity problem falls into the scope of the results of Kondrat'ev's "parabolic" paper 
[4~8] . The case of "expanding paraboloids" as t — > 0~, similar to the above cases for the 
bi-harmonic equations and others, was not treated before by the same reasons (the lack 
of proper spectral theory and matching via a boundary layer approach, which are very 
difficult to fully justify). 

Here, according to the principle (11.81) . we very briefly discuss a new types of operator 
pencils that is needed to tackle other regularity problems. The simplest such model 
satisfying (11. 8p is the linear fourth- order hyperbolic (wave) equation 

(8.44) u tt = ~u xxxx in Q , 

also known as the ID linear beam equation. We consider the same Dirichlet problem 
as in (14. ip with two bounded initial functions u(x, — 1) = u$(x) and u t (x, — 1) = ui(x). 
Concerning the fundamental solution and necessary spectral properties, we follow [29|. 
Thus, the fundamental solution of (I8.44p has the self-similar form 

(8.45) b (x,t) = Vt F(y), y = f t , where b {x,0) = 0, b ot (x, 0) = S(x) 

in the sense of bounded measures. The rescaled kernel Fq = Fo(\y\) is symmetric and 
solves the ODE 

(8.46) BF = -F^ - \ F"y 2 - \ F'y + \ F = in R, J F = 1. 
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Integrating (18.461) once yields 

oo 

(8 .4 7) _ F >»_i F > y 2 + i F y = o F{y) = _L J mi^M dz . 



A WKBJ-type asymptotic analysis of the ODE ( 18.471) yields the behaviour like ( 11. 81) : 

(8.48) F(y)~C 1 y-#cos(^ + C 2 ) as y -> +oo, d ^ 0. 

Similar to (14. lip , the first scaling is 

= v(y,r), y=-rs, r = -ln(-t), 
(8-49) ^ | 

W TT + W r + = = -Vyyyy ~ J y^Vyy ~ J yVy. 

Not paying enough attention to the functional setting of the operators involved in weighted 
L 2 -spaces |29j . we concentrate on the polynomial eigenfunctions of a quadratic pencil to 
appear; see Markus [63] for necessary concepts and theory of linear operator pencils. 
Namely, to find eigenfunctions, we set 

(8.50) w(y,r) = e^r k (y) C*(A fc )^ = B*P h - + h)Pk ~ A fc (^)' V = 0. 
Looking for finite polynomial eigenfunctions (recall that one needs two sets of eigenfunc- 
tions ty* = {ipk(y)i ^tiv)} as we U as the adjoint ones \& = {ipk(y), 4>k(y)} [29])- 

(8.51) i } k{y)^ ^kiv) = y k + ••• are ^th-order polynomials, 

substituting into (18.501) and keeping the higher-degree terms ~ y k yields the following 
quadratic equation for eigenvalues: 

(8.52) X 2 k + (k + l)X k + + f = 0. 
This gives two, shifted by -1, series of eigenvalues, 

(8.53) K = -\ and K = ~l~ l for k = 0,1,2,.... 

Let us specify a half of those generalized Hermite polynomials, which is a complete closed 
set. There is another one obtained in a similar manner, which is obviously necessary for 
the second-order in t PDE ( I8.44p . 

Proposition 8.1. The eigenfunctions of the pencil with the spectrum o~ + = {X~^} 

in <\8.53§ are given by the normalized polynomials 

(8-54) r k (y) = ^[y k + J2%^(y k ) m ], fc = 0,1,2,.... 

Assuming as usual that <p(r) is a slow growing function (so (14.71) . (14.81) . ... hold), we 
introduce the variable z in (17. ip and the resulting BL- variables as in (I7.2p . where 

(8.55) i = ^{r){l-z) 1 z = 1 fc, D y = -<p(r)D $ . 

Assuming that the class of generic solutions admitting a BL-representation is under 
scrutiny and performing the cut-off as in (I7.14p . where v TT = v TT H, by the expansion 
(17.201) (we currently take into account the first coefficient only, with (1,F) = 1 by the 
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orthonormality, over the first half (I8.54p of all the polynomials, and this defines the class 
of solutions involved), so eventually, by the asymptotics (18.481) yields 

(8.56) < + o(, + ... ~ 7 iC 3 cwi cos(^ + C 4 ) +... (§ = 3-i|), 

where we omit all smaller terms containing ip', if", and others. For slow varying coefficients 
a (r), the leading term is indeed a' Q , so, in this critical case, the governing ODE is 

(8.57) J ~ 7 C>§ cos(^ + C 4 ) + ... In \a {r)\ ~ j cos + C 4 ) dr. 

Clearly, unlike ( 18.211) . by setting <p 2 {t) = s, we see from ( 14.7ft that a converging integral for 
the irregularity is impossible. However, this does not mean a "total" regularity. Actually, 
this again implies that the regularity issue here is rather subtle and, as usual for oscillatory 
kernels, can be achieved by a proper oscillatory cut-off of suitable slow growing f(r). For 
general ip(r), a more delicate analysis of the ODE (18.561) with all the terms included is 
necessary, which leads to some technical questions. A full justification of the general 
approach then becomes difficult and even questionable. 
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Appendix A. Comment: On using majorizing order-preserving operators 



It is clear that the main difficulty in establishing optimal regularity/irregularity criteria for 
higher-order linear parabolic PDEs such as (|4.ip . (jl,4p (or third-order dispersion ones in (|8.5p and 
(I8.25h to be studied later on) is the fact that these flows do not obey the Maximum Principle and 
do not exhibit order-preserving features. As is well-known, the latter ones have been key for the 
heat equation (|3.ip and other nonlinear second-order diffusion-reaction- absorption-convection-... 
equations; see references in [T| |2| 1361 176] , 

Therefore, before going into eigenfunction expansion techniques for solving the regularity 
problem for (14. ip , we cannot avoid a temptation to briefly discussing a rather formal opportunity 
to deal with the traditional techniques based on the Maximum Principle even for the higher- 
order parabolic flow. This is about the idea of majorizing order-preserving flows for higher-order 
parabolic equations [31j . which we explain relative to the poly-harmonic flow (11. 4p . Namely, 
given the rescaled oscillatory kernel (|5.5p . we construct the corresponding majorizing kernel 
F(y) > such that 

(A.l) \F(y)\ < D m F(y) in R N , where JF = 1, 

and D m > 1 is a constant called the order deficiency (or order defect) of this majorizing kernel. 
It is clear that the optimal order deficiency is {D\ = 1 since F > 0) 

(A. 2) = f \F\ > 1, corresponding to the non-smooth kernel F(y) = -r4- \F(y)\ 

m 

(sometimes, dealing with not C°° kernels F{y) is not convenient, especially when determining 
its eigenfunctions by differentiation via (|5,16p . with F \-t F, etc.). 

We then introduce the corresponding majorizing order-preserving integral evolution equation 
in 1H N x R + (it does not have a PDE representation in general) defined by the formal convolution 

(A.3) u(t) = b(t) * u , with u (x) > 0, where b(x,t) = t~^F(y), y = -^. 

The corresponding comparison theorem is as follows |31| : 

(A.4) \u (x)\ < jj^ u (x) in R N => \u(x,t)\ < u(x,t) in R N x E + . 

The proof is straightforward in view of (jA.ip via comparison of the solution of (|A.3p with the 
parabolic flow 

(A.5) u(x,t)=b(x-;t)*u (-) = t-^ jF(f^)u (y)dy in R N x R + . 

Indeed, subtracting (|A.3|) and (|A.5|) yields for the difference the desired inequality: 
w(x, t) = u(x, t) — u(x, t) = t~ 2^r f \F(-) uq — F(-) uo] 

(A.6) 

> / [F(-) n - ^ |n-)l So] = J [ DmF{ t im ^ 0. 

by the definition of the majorizing kernel (|A.1|) . Since F > 0, i.e., the integral majorizing 
evolution (IA.3P is order-preserving, solutions u(x, t) admit constructing super-solutions U(x, t) 
in the standard sense, i.e., the following comparison holds: 



(A.7) U{t)>b(t)*u => < u(x, t) < U(x, t) in 
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In [31], such comparison ideas resolved some asymptotic problems for semilinear parabolic equa- 
tions related to Fujita critical exponents and around. 

It is clear that to construct a majorizing comparison theory for the IB VPs in domains such as 
Qo, further estimates of "Green's function" corresponding to the "majorizing non-translational 
flow" (|A.3P hi Qq (as usual, such an estimate is expected to be related to (|A.1[) ) are neces- 
sary. These are difficult questions, which are not discussed here. We just mention that such 
a true possibility would allow, in a reasonable simpler (and rigorous!) way, to get a proper 
upper estimates on solutions in Qq to prove that u(0, 0) = 0, which would guarantee its regular- 
ity, by constructing proper super-solutions. As a formal clue, it would be very efficiently and 
attractively easy to use Petrovskii's ansatz (13.5|) . 



2m 



where we have replaced the Gaussian structure in (|3.5I) by that in the estimate ([5.7]) . which 
is supposed to be true for the majorizing positive kernel F{y). As we will show (cf. (|7,33p 
for m = 2), the space-time structure (|A.8j) would give a correct critical behaviour for R(t) as 
t —* 0~ , but we do not know if this can be justified in such a comparison way. (an open problem, 
but an extremely attractive one). We must admit that, besides some other fundamental open 
questions, in any case, unlike (13. 51) . using (1A.8|) for any m > 2 assumes complicated integral 
calculus to establish the validity of the integral inequality in (|A.7p . More importantly, such an 
integral majorizing approach cannot be used for any lower bound on u(x, t) to get the irregularity 
counterpart, when \u(0, 0)| > 0. As one can see, corresponding alternative "minorizing" theory 
in similar lines cannot be constructed in principle. Therefore, we had no intention and in 
fact necessity to check whether functions such as (IA.8P can satis fy ( I A. 71) or similar integral 
inequalities. However, we state the following principal open problem 21 ]: Is there any hope to get 
sharp regularity results by using majorizing-like calculus associated with some integral operator 
in Qq with a positive kernel? 
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Personally, the author suspects that this is not possible, especially since the "cut-off" must be then 
inherited by the method (how and why?); though then proving the negative result (a principal non- 
applicability of such barrier techniques for m > 2) would be also important. 
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